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Introduction 


Nous allons, cette année, aborder dans le cours de Mécanique du solide et 
Mécanique analytique la formulation moderne des principes de la mécanique des 
systèmes dynamiques à un nombre fini de degrés de liberté. Le Lecteur consultera 
avec profit les ouvrages classiques [1, 10] qui ont inspiré ce cours. 

Le formalisme mathématique de la mécanique rationnelle développé par Joseph- 
Louis Lagrange (1736-1813) dans un corpus scientifique considérable, notamment 
sa Mécanique analytique (1788), a conduit à une généralisation des principes de 
la mécanique newtonienne à des systèmes dynamiques plus élaborés que celui du 
simple point matériel. Le formalisme lagrangien que nous allons introduire tire son 
origine du principe des travaux virtuels (d’Alembert, Maupertuis) qui a recours 
à la notion de mouvements virtuels d’un système pour déterminer le mouvement 
. . . réel de ce système. Dans le cas de systèmes soumis à des forces conservatives, 
les équations de Lagrange (équations du mouvement) sont dérivées d’une seule et 
unique fonction, le lagrangien L, sans avoir à prendre en considération les forces 
de liaisons (holonomes) souvent très complexes. D’où une simplification concep- 
tuelle et pratique de la mise en équation des problèmes mécaniques. Ce formalisme 
est également géométrique car indépendant du choix d’un système de coordonnées 
(généralisées) ; d’où une extension naturelle au cas d’espaces de configuration très 
généraux (variétés différentiables). Nous émaillerons cette partie du cours de nom- 
breux exemples illustratifs, notamment le problème des N corps, les pendules, cer- 
tains systèmes à liaisons holonomes, le couplage minimal d’une particule chargée à 
un champ électromagnétique extérieur, etc. 
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L’autre approche de la mécanique des systèmes que nous aborderons concerne le 
formalisme hamiltonien (Sir William Rowan Hamilton, 1805-1865) introduit dans 
une série de travaux, notamment dans son article “On a General Method in Dyna- 
mics” (1834). Alors que le formalisme lagrangien mettait en jeu une fonction L de 
l’espace tangent à l’espace de configuration du système (espace des couples position- 
vitesse), le formalisme hamiltonien a également recours à une unique fonction, H, 
mais définie cette fois-ci sur V espace cotangent à l’espace de configuration (espace 
des couples position-impulsion) ; cette fonction est Y hamiltonien du système et cor- 
respond à Y énergie ou générateur de l’évolution temporelle du système. Là encore, 
la théorie hamiltonienne admet des généralisations géométriques (variétés symplec- 
tiques, variétés de Poisson). Les méthodes modernes de quantification (décrivant le 
passage d’une description classique à une description quantique de l’univers physi- 
que) utilisent le formalisme hamiltonien de manière fondamentale. 

Ces différents aspects de la mécanique analytique trouvent naturellement un 
champ d’application dans le chapitre important de la mécanique rationnelle que 
constitue la mécanique du solide (Euler, Poinsot, Lagrange, Kovalewski, etc.). Les 
mouvements du corps rigide (par exemple une toupie) sont très riches et leur étude 
subtile. Le fait que le corps solide ne présente pas, en général, de symétries parti- 
culières conduit à la notion importante <T opérateur d’inertie servant à décrire ses 
différents mouvements en présence ou non de forces extérieures. La configuration 
d’un solide sera, nous le verrons, déterminée par les éléments d’un groupe, le groupe 
euclidien SE(3), de l’espace euclidien tridimensionnel : la position d’un point origine 
du solide et l’orientation générale du solide par rapport à ce point. L’étude des chan- 
gements de référentiels (passage Laboratoire-Solide) paramétrés par le temps est une 
étude obligée et riche d’enseignements (mécanique dans les systèmes non inertiels ; 
par exemple l’étude du pendule de Foucault, des ouragans et typhons, etc.). La dy- 
namique d’un solide libre avec point fixe sera étudiée grâce aux théorèmes généraux 
de la mécanique et aussi dans le cadre hamiltonien (équations d’Euler). Le problème 
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de la toupie avec un point fixe, plongée dans le champ de pesanteur, est associé au 
nom de . . . Lagrange : l’étude de certains des mouvements de la toupie de Lagrange 
sera également abordée. 

Le champ d’application de l’étude du corps solide est vaste et important non 
seulement sur le plan de la mécanique abstraite 1 (alias analytique) mais aussi et 
surtout en mécanique appliquée où sont à l’œuvre les phénomènes gyro scopiques. 2 
Citons, à titre d’exemples, les gyrocompas ou boussoles gyroscopiques déterminant 
le nord (deux degrés de liberté), les gyroscopes servant à stabiliser les satellites, à 
déterminer l’horizon artificiel dans les avions, etc. La stabilisation des (moto) cyclistes 
résulte aussi de l’effet gyroscopique. Rappelons enfin que l’explication du phénomène 
de précession des équinoxes (précession de l’axe de rotation de la terre — toupie apla- 
tie sous l’effet des forces de marées dues au soleil et à la lune — avec une période 
de 25800 ans autour de la perpendiculaire au plan de l’écliptique) relève encore de 
la mécanique du solide. 

Mes plus sincères remerciements vont à Jean-Philippe Michel pour sa lecture 
attentive et critique du manuscrit de ce cours. 


1. Certaines toupies constituent des exemples de systèmes dynamiques intégrables, i.e. résolubles 
par quadratures, et un champ de recherche privilégié en mathématiques et physique mathématique ; 
l’aspect quantique de ces systèmes intégrables est un objet d’étude actuellement très actif. 

2. En l’absence de moment de forces extérieures, un solide possède un moment angulaire 
constant ; d’où l’importance des dispositifs mettant en jeu des gyroscopes pour la stabilisation 
des véhicules terrestres, maritimes et aériens. 
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Chapitre 1 


Les équations de Lagrange 

1.1 Une introduction heuristique 

Considérons pour (bien) débuter les trajectoires r (t) possibles d’une particule de 
masse m se déplaçant sous l’influence d’un champ de forces extérieur F dérivant, par 
exemple, d’un potentiel V (r) indépendant du temps t. Les équations du mouvement 
de Newton s’écrivent, on le sait, 

mr{t) = F(r(f)) avec F = — grad(F) (1.1.1) 

où grad(F) = dV/dr désigne le gradient 1 de la fonction différentiable V. 

Désignons par T = |m||r|| 2 l’énergie cinétique de la particule; elle varie, bien 
sûr, au cours du temps le long de chaque trajectoire r(f). Le potentiel étant ici 
stationnaire, i.e. dV/dt, l’énergie totale est une constante du mouvement : 

H — T + V — const. (1.1.2) 

Exercice 1.1.1. Vérifier la loi de conservation de l’énergie (1.1.2). 


1. Abus usuel de notation : on devrait écrire 


grad(F) = — 

où la barre désigne la transposition; on note vw = (v, w) le produit scalaire ordinaire (euclidien) 
de v, w G K 3 et ||v|| = \J (v, v) la norme de v G K 3 . 


1 
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Considérons, maintenant, la nouvelle expression 

(1.1.3) 

définie par l’étrange différence de l’énergie cinétique, T, et de l’énergie potentielle, V . 
Cette expression est clairement une fonction L( r, r) de 1a. position r et de 1a. vitesse r 
de la. particule. On trouve facilement 


L = T — V 


d L 

— — = mr et 

or 


dL _ _dV 

dr dr 


de sorte que les équations de Newton (3.2.22) peuvent se réécrire de la. manière 
suivante 

mr 


dV 


dL 


7 ... 

or 


dr dt \ dr J 
Nous avons donc prouvé le résultat suivant : 


Définition-Théorème 1.1.1. Soit T l’énergie cinétique d’une particule plongée 
dans un potentiel V ; on appelle lagrangien du système la fonction L = T — V. 
Le système des équations du mouvement de Newton est équivalent au système des 

équations de Lagrange 

( 1 . 1 . 

Remarque 1.1.2. Le théorème précédent reste, bien entendu, valable dans le cas 
général d’un potentiel V(r,t) dépendant explicitement du temps. 

Exercice 1.1.3. Soit L( r,r) = Om(||r || 2 — ca 2 ||r|| 2 ) un lagrangien défini sur M 3 x M 3 ; 
écrire les équations de Lagrange. En donner la solution générale. Interprétation 
physique ? 

Exercice 1.1.4. Soit L(r,9,r,9 ) = |(r 2 + r 2 9 2 ) ; écrire les équations de Lagrange. 
Donner la solution générale ( r(t ), 9{t)) de ce système d’équations différentielles. Que 
représentent ces courbes du plan euclidien ? 

Exercice 1.1.5. Trouver l’expression du lagrangien L{9,9) d’un pendule simple de 
masse m, longueur I dans le champ de pesanteur g = const. 


d , 

(dL\ 

dL n 

& 

dr 

dt 

\ Or J 

~ 1k ~ 

—r~ — r 

dt 
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1.2 Illustration : Equations de Fermât 

Un exemple important utilisant les équations de Lagrange concerne l’optique 
géométrique dans la formulation qu’en a donnée Pierre de Fermât (1601-1665). Selon 
le Principe de Fermât, les rayons lumineux se propagent dans un milieu d’indice 
de réfraction n(r) variable selon des courbes minimisant le chemin optique entre deux 
points quelconques. Ce principe revient, on le verra, à décrire les rayons lumineux 
par les solutions des équations de Lagrange pour le lagrangien 

L(r, r) = n(r) ||r ||. (1.2.1) 


Ce lagrangien s’interprète comme le chemin optique élémentaire , ns, parcouru par 
la lumière dans un milieu d’indice n à la vitesse (scalaire) s = ds/dt ; ici t est un 
paramètre décrivant les rayons lumineux (courbes de l’espace euclidien E 3 ) et s 
désigne l’abscisse curviligne définie par la primitive suivante 


s = 



dt. 


Nous noterons 


la “vitesse unitaire” . 2 


d.r r 

ds || r 


Les équations de Lagrange (1.1.4) prennent alors la forme 


d /<9L\ 
dt \dr J 


d L 

dr 




0 


ou encore 


d ds 

— nu — 
ds dt 


grad n 


ds 

dt 


puisque ds/dt = ||r|| > 0 (on ne considère pas les points de rebroussement). 


2. Attention : la vitesse de propagation de la lumière est bien infinie ( !) dans le cadre de 
l’optique géométrique si t désigne le temps galiléen. Dans cet exemple t est, soulignons-le, un 
paramètre arbitraire servant à décrire les rayons lumineux. Rien à voir avec le temps . . . 
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Proposition 1.2.1. Les équations de Fermât gouvernant l’optique géométrique sont 
données par les équations de Lagrange pour le lagrangien (1.2.1) et prennent la forme 

( 1 . 2 . 2 ) 

Corollaire 1.2.2. Dans le vide, n = 1, les rayons lumineux empruntent les droites 
(géodésiques) euclidiennes d’équation paramétrique 

r(s) = r(0) s + r(0) 


d 


dr 

— nul = gradn 

ds 

& 

d^ = U 


pour tout sel. 


Exercice 1.2.3. Supposons que l’indice de réfraction soit, dans le plan, donné par 
la fonction discontinue 

n ( x v ) = ! ni (y > °) 

[ ,V) \ n 2 (y< 0) 

où ni et n 2 sont deux constantes (positives). Déduire des équations (1.2.2) les Lois 
de Descartes 3 

ni sin R = n 2 sin i 2 ( réfraction ) 
ii = —ï 2 ( réflexion ) 

où ii et i 2 sont les angles orientés formés par les rayons lumineux et la normale e y 
au dioptre y = 0. A f.B. Ne pas chercher à déterminer le gradient de l’indice de 
réfraction (car cet indice n’est pas une fonction continue!). 


Exercice 1.2.4. Déterminer les trajectoires des rayons lumineux dans le demi-plan 
supérieur, H + = {( x,y ) G M 2 | y > 0}, si l’indice de réfraction est n(x,y ) = l/y. 


1.3 Equations de Lagrange 

1.3.1 Formalisme intrinsèque 

Nous avons introduit les équations de Lagrange (1.1.4) dans un cas réellement 
très particulier (cas d’une particule non relativiste dans un potentiel extérieur) et, de 
plus, dans un système de coordonnées spécial (les coordonnées cartésiennes de M n , 
pour n = 1, 2, 3, . . .). Les questions suivantes viennent alors naturellement à l’esprit : 


3. Ces lois sont souvent attribuées à Snell & Descartes. 
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1. Peut-on généraliser les équations de Lagrange au cas des systèmes quelconques 
à un nombre fini de degrés de liberté (e.g. le problème des N corps en inter- 
actions mutuelles) ? 

2. Les équations de Lagrange possèdent-elles un caractère intrinsèque (indépen- 
dant du système de coordonnées choisi sur l’espace de configuration) ? 

La réponse à la première question sera apportée par les nombreux exemples qui 
émaillcront la suite de l’exposé. Quant à la deuxième question, la réponse est fournie 
par la 

Proposition 1.3.1. Soit (q 1 , . . . , q n ) un système de coordonnées arbitraire sur l’es- 
pace de configuration d’un système mécanique. Les mouvements de ce système sont 
donnés par les solutions des équations de Lagrange 

(1.3. 

avec L = T — V , où T désigne l’énergie cinétique et V l’énergie potentielle du 
système. 


d 

( d _L\ 

ÔL 

dt 

\dq l ) 

dq i 


défi 

dt 


Mi — 1, . . . , n 


Démonstration. Nous nous limiterons à un système à un degré de liberté, n — 1. La 
généralisation au cas n > 1 est laissée en exercice. 

Considérons donc un changement de coordonnées arbitraire q i— > q* = Q(q) où 
Q : M — > M est une application monotone (un difféomorphisme local), i.e. vérifiant 
partout dq* /dq = Q'(q ) fi 0. Le lagrangien s’exprime ainsi dans chacun des systèmes 
de coordonnées selon L = f(q,q) = f*(q*,q*)- Notons que pour toute courbe q(t) 
on a dq*/dt = d(Q(q))/dt = Q'{q) dq/dt — dérivée d’une fonction composée. La 
transformation des coordonnées de configuration et de vitesse est, en définitive, 



i— ► 



= Q(q) \ 
= Q'{q)q ) 


(1.3.2) 
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On trouve aisément, cf. (1.3.2), 

dL 

dq 


dL dq* dL dq* 

dq* dq dq* dq 


puisque dq* /dq = 0. D’autre part 

dL 
dq 


dL dq* dL dq* 
dq* dq dq* dq 

dï Q{q) + W Q(q)q 


On obtient donc 


d_ 

dt 


dL 

dq 


dL 

dq 


Jt (w) Q,(q) + W Q " (q) q -W Q ' (q) -W> Q " (q) q 
Q\q) 


' d 

(dL\ 

dL] 

dt 

\dq* J 

1 * 
|c^> 

1 


Si les équations de Lagrange (1.3.1) sont vérifiées dans les coordonnées (q, q), elles 
le seront automatiquement dans tout autre système ( q*,q *) obtenu par un difïéo- 
morphisme (1.3.2) découlant de la condition Q\q) ^0. □ 


1.3.2 Exercices illustratifs 

Le but de ce sous-chapitre est de fournir, sous forme d’exercices relativement 
détaillés, des exemples de problèmes physiques mis en équation grâce au formalisme 
lagrangien. 

Le système newtonien des trois corps 

Considérons le système de trois corps Mi, M 2 et M 3 , de masses mi,ni 2 et rris, en 
interaction gravitationnelle mutuelle. Le potentiel d’interaction newtonienne entre 
les points et M } (avec i É j) es t donné par 

Grnjrij 


bp'( r i, r 3 , r 3 ) 


r * - r i 
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où r, : = OMj représente la position de l’astre M t (i = 1, 2, 3) par rapport à une 
origine arbitraire O, et G la constante de Newton. On désigne par r j la vitesse du 
point Mj à un instant donné. Justifier le lagrangien du système 


£(ri,r 2 ,r 3 ,ri,r 2 ,r 3 ) = ^mi||iq || 2 + ^m 2 ||r 2 || 2 + ^m 3 ||r 3 || 2 


Gm 2 m 3 Gm 3 mi Gmim 2 
"fin 1 1 J - 1 1 1 1 J~ 


r 2 - r 3 


r 3 - r 3 


ri - r 2 


en vérifiant que les équations de Lagrange sont bien équivalentes aux équations du 
mouvement de Newton pour ce système. 


Le système général des N corps 

Nous envisageons maintenant le cas général d’un système de points matériels 
Mi, . . . , Mtv de masses mi, . . . ,rriN en interaction mutuelle et soumis à des forces 
extérieures. Nous supposerons que toutes les forces enjeu sont conservatives. Désignons 
par Vij le potentiel d’interaction entre Mi et Mj (avec i ^ j ) et supposons qu’il ne 
dépende que de = || r.j — r ^ 1 1 , la distance relative de Mi et Mj. On notera Vj ext le 
potentiel dont dérive la force extérieure F^ xt appliquée au point Mj. 

Justifier que 

= Fy fe - &] 


pour tout k = 1, 2, 3, où 


F -J-Yti 

ij ~ dn 


représente la force à laquelle est soumise Mi interagissant avec Mj. Vérifier que la 
troisième loi de Newton est vérifiée si Vj 3 = Vji. 

Considérons le lagrangien suivant 


n n 1 N N 

L ( r >r) = - J]mj||rj|| 2 - - V v 


ext 


(1.3.3) 


i— 1 


i. 7=1 


i= 1 


où r = (ri,r 2 ,r 3 ) et r = (fi,r 2 ,r 3 ). Montrer que les équations de Lagrange (1.3.1) 
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restituent les équations du mouvement du système des N corps en interaction, 

fi = V F, + Ff*. 


pour tout i — 1, . . . , N. 

Le pendule double 

Ce système est constitué de deux pendules M\ & M 2 , de longueurs i\ & 1 2 
constantes et de masses rrii & m 2 plongés dans le champ de gravitation g = ge x 
avec g = const. > 0 ; le point de suspension du deuxième pendule est le point M\ 
et celui du premier, l’origine O fixe du système de coordonnées cartésiennes (x, y) 
dans le plan euclidien. On désigne par 9\ (resp. 9 2 ) l’angle que forme le pendule M\ 
(resp. M 2 ) avec la verticale. 

Déterminons le lagrangien du système. Posons M n = O + ( x n ,y n ) pour n — 1,2 
dans le repère orthonormé direct (e x e y ) du plan pendulaire ; on a alors X\ = i\ cos 9 1 , 
yi = I\ sin 6i & x 2 = I\ cos 6 ± + I 2 cosd 2 , y 2 = I\ sindi + I 2 sind 2 . 

L’énergie cinétique de AR est T x = ^m\{x\ + ÿf) = celle de M 2 est 

alors T -2 = \m 2 (x \ + ÿ|) = \m 2 (I\dl + I\Q\ + 2Q£ 2 cos(di - 9 2 )9i9 2 ). 

Les énergies potentielles sont, de même, Vf = —migxi = —m\g£\ cos 9\ ainsi que 
V 2 = —m 2 gx 2 = —m 2 g{£\ cos9\ + £ 2 cos 9 2 ) — à des constantes additives près. 

Le lagrangien total est alors L — T\ + T 2 — V\ — V 2 , c’est-à-dire 

L(9i, 9 2 , 9\, Ù 2 ) = -(mi + nr 2 )£^d^ + —m 2 £. \9 2 + rn 2 lil 2 cos(di — 9 2 )9i9 2 
+(?7ii + m 2 )gli cos 9\ + m 2 g£ 2 cos 9 2 . 

Les équations de Lagrange (1.3.1) fournissant les équations du mouvement du 
système s’écrivent ici d(dL / 39 n ) / dt — dL/d9 n = 0 pour n = 1,2. On obtient ainsi, 
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après quelques simplifications, le système suivant 

0 = (m 1 + m 2 )£ îdi + m 2 £ 2 cos(0i - 0 2 )0 2 + m 2 t 2 sin (0 X - 9Q9\ (1.3.4) 

+(mi + m 2 )g sin 9 1 

0 = £ 2 9 2 + £\ cos(0i — 9 2 )9i — I\ 8111(0! — 0 2 )$i (1.3.5) 

+g sin0 2 . 

Ce système d’équations différentielles non linéaires couplées n’est pas intégrable 
analytiquement. On sait par contre, depuis Henri Poincaré, que ses solutions ex- 
hibent un comportement chaotique , c’est-à-dire, une sensibilité structurelle aux condi- 
tions initiales (impossibilité de retrouver les mêmes trajectoires du système en 
répétant l’expérience — avec des conditions initiales identiques sur le plan expérimental). 

Nous pouvons, par contre, étudier les petits mouvements du système autour de 
la position d’équilibre (stable!) évidente 4 

0i = 0 2 = 0. 

On trouve alors aisément le système linéaire d’équations différentielles du second 
ordre couplées suivant 

0 = 9\ + ji—9 2 —9i (1.3.6) 

ii 1 1 

0 = 0i + — 0 2 + ^-0 2 (1.3.7) 

où 

m 2 

9 = — : — 

mi + m 2 

est tel que 0 < p < 1. 

Exhibons les fréquences propres naturelles du système dans le cas simple où 

h = £ 2 = £ 


4. On vérifie que 9\{t) = 9 2 {t) = 0 est une solution exacte des équations (1.3.4) et (1.3.5). 
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en recherchant des solutions particulières du système (1.3.6) et (1.3.7) sous la forme 


(t) = A cos eut 
d 2 (t) = B cos ut 

où A, B et lv > 0 sont des constantes à déterminer. On trouve aisément 


u 


u 


B 

. 1+Al i 

B 

1 + I 


9_ 

I 

g_B 

IA 


(1.3.8) 


et, en divisant membre à membre les deux équations précédentes, (B /A) 2 = 1/p. 
En reportant maintenant ce résultat : 


± VU 

dans une des deux équations (1.3.8), il vient alors 

2 = 9/t 

± i±VU' 

Les deux pendules oscillent en phase aux basses fréquences 


9f(t) = Acos(uq_f) 
0t ( t ) = ~^= cos (u+t) 

y/U 


& 


= 


g/t 


1 + sJ Ji 


ou en opposition de phase aux hautes fréquences 


0 l ( t ) = AcosQjH) 

A 

9 2 if) — cos(cu_t + 7r) 

VU 

où 9f (0) = A reste arbitraire mais . . . petit. 


& 


OJ _ = 


g/t 


i - VU 


(1.3.9) 


(1.3.10) 


1.3.3 Liaisons holonomes 

Si un point matériel est astreint à se mouvoir sur une courbe ou une surface 
dans l’espace euclidien (par exemple sur une sphère de rayon donné dans le cas du 
pendule sphérique), on dit que le système présente une liaison holonome. 
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Donnons maintenant la définition générale de la notion de liaison holonome. 
Considérons un système mécanique à n degrés de liberté, i.e. dont la configuration est 
déterminée par n paramètres indépendants. Par exemple un système de N particules 
évoluant dans l’espace euclidien E 3 est un système à n — 3N degrés de liberté ; dans 
ce cas (E 3 ) N = E 3 x ... x E 3 est l’espace de configuration. 5 Supposons maintenant 
que les points soient contraints à se déplacer sur une hypersurface définie par une 
équation 

F(Mi, . . . , M n ) = 0 (1.3.11) 

où F : (E 3 ) n — » W 71 est une application différentiable (avec 0 < m < n) définissant 
la contrainte ; nous définirons le nouvel espace de configuration par 

Q = F- 1 ( o). 

Cette hypersurface (ou “variété” ) est décrite par n — m coordonnés “généralisées” ; 
dans ce cas dirn(Q) = n—m est le nombre de degrés de liberté du système contraint. 6 
Nous dirons que la condition (1.3.11) est une liaison holonome pour le système. 

Remarque 1.3.2. Les liaisons holonomes ne mettent pas enjeu des conditions sur 
les vitesses des points du système. Par exemple la liaison décrivant une roue roulant 
sans glisser sur une route est une liaison non holoriome — elle exprime le fait que 
la vitesse du point de contact de la roue avec le sol est nulle. 

Exemple 1.3.3. Un haltère formé de deux masses reliées par un manche de longueur 
I > 0 a pour espace de configuration 

Q = {(ri, r 2 ) G M 3 x M 3 | F(r 1; r 2 ) = \\r 1 - r 2 \\ 2 - f = 0} 

et cet espace de dimension 5 a la topologie Q = M 3 x S 2 où S 2 désigne la 2 -sphère 
(de rayon £ > 0) plongée dans M 3 . 

5. Il faudrait, en fait, considérer que l’espace de configuration est plutôt ( E 3 ) N — {collisions}. 

6. Ceci découle du fait que F est supposé de rang maximum (un point technique important 
mais délicat que nous ne développerons pas). 
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Exercice 1.3.4. Donner la fonction F : M 2 x M 2 — > M 2 définissant l’espace de confi- 
guration Q du double pendule plan. Vérifier que Q = S 1 x S 1 ( tore 2- dimensionnel) . 

Comment le formalisme lagrangien se présente-t-il dans le cas de systèmes avec 
liaisons holonomes ? Nous allons montrer que les équations de Lagrange se formulent 
de fait en oubliant complètement les forces de liaison (contrairement au formalisme 
newtonien, plus difficile à mettre en œuvre dans un sens). 

Illustrons ceci en donnant un exemple élémentaire où un point matériel est as- 
treint à se déplacer sur une courbe du plan euclidien (par exemple un cercle de rayon 
donné pour le pendule simple). Nous supposons, en fait, notre point soumis à un 
potentiel de confinement (potentiel harmonique intense dans la direction orthogonale 
à la courbe créant une force de rappel intense vers les points de la courbe). 

Soit r = (x, y) G M 2 la position d’un point matériel de masse m dans le plan (on 
a choisi un système de coordonnées euclidiennes) et soit 

V u (x, y) = V 0 (x, y) + ^ mu 2 y 2 (1.3.12) 

le potentiel tout à fait général dans lequel est plongé ce point de masse m. Nous 
allons ensuite considérer la limite u — > oc qui permettra de confiner le point matériel 
sur l’axe des x. 

Le lagrangien du système est alors L(x, y , x, ÿ) = fim(x 2 +y 2 ) — V U! (x, y). L’équation 
de Lagrange pour y s’écrit d/dt(dL/dy) — dL/dy = m(ÿ + u> 2 y) +dVo/dy = 0. Dans 
la limite des grandes pulsations u le dernier terme sera tout à fait négligeable et 

ÿ + u 2 y 2 « 0. 

Cette équation différentielle possède, on le sait, une intégrale première (reliée à 
l’énergie) e = ÿ 2 + u 2 y 2 = w 2 r/^ ax = const. Donc -y max < y < y max où y max = y/ë/t o. 
Finalement y max — > 0 si uj — ■» oo et le point est donc confiné sur la courbe y = 0 par 
cet artifice mathématique. 

L’équation de Lagrange pour x s’écrit dj dt(dL / dx) — dL / dx = mx+dV 0 /dx = 0. 
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Dans la limite ca — > oo nous aurons finalement 


mx + V'(x) = 0 où V(x) = Vq{x, 0) (1.3.13) 


puisque lim^oo y — 0. Nous venons de prouver le résultat suivant 

Proposition 1.3.5 ([1]). Un point matériel est soumis au potentiel (1.3.12) de 
confinement V ÙJ (x,y) dans le plan euclidien. Soit (x u (t), y w {t)) la solution générale 
des équations de Lagrange avec conditions initiales (^(0) = 0, ÿ w (0) = 0). Lorsque 
u — » oo il existe une limite x(t) = lim^oo x u (t) qui vérifie l’équation de Lagrange 


d / d L q 
dt V dx 


dL 0 

dx 


= 0 


ou 


U() L\y = y = Çj. 


Ce résultat se généralise immédiatement au cas de toute courbe dans le plan et, 
plus généralement, au cas de toute hypersurface d’un espace de configuration. 

En fait, les équations de Lagrange pour un lagrangien L d’un système défini par 
des liaisons holonomes (1.3.11) se ramènent aux équations de Lagrange 
d ( d L 0 \ dL 0 


, . r... , „ . = 0 où Lq — L |p_p_ n (1.3.14) 

dt \ défi ) dq 1 \F-F-o \ J 

pour i = 1, . . . , n — m exprimées dans un système de coordonnées (q 1 , . . . , g n_TO ) de 

la surface Q = F^ 1 (0). 


Exercice 1.3.6. Ecrire les équations du mouvement d’un pendule sphérique de lon- 
gueur é dans le champ de pesanteur g = const. N. B. On utilisera, bien sûr, un 
système de coordonnées sphériques ( 6,ip ) de S 2 . 

Démonstration. La position d’un point sur la sphère S 2 C M 3 de rayon I est para- 
métrée par r = (x, y , z) = £(sin 6 cos <p, sin 6 sin cos 9), en coordonnées sphériques ; 
ici p> G (0, 27t) désigne la longitude et d G (0,7 r) la colatitude. La vitesse du point 
mobile sur la sphère étant v = r, il vient aisément ||v|| 2 = £ 2 (d 2 + sin 2 9fi 2 ). Comme 
l’énergie potentielle est V = mgz = mg£cos0, à une constante additive près, on ob- 
tient le lagrangien L = |m||v|| 2 — V du pendule sphérique en terme des coordonnées 
sphériques et des vitesses associées 

L(O,(p,0,(p) = ^m£ 2 (9 2 + sin 2 9 fi 2 ) mgicosd. 
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Les équations de Lagrange se lisent maintenant 

9 — sin 9 cos 9 p 2 — ^ sind = 0 (1.3.15) 

[sin 2 9(p\ = 0 (1.3.16) 

LLL 

ce qui implique une constante du mouvement (reliée au moment angulaire ) 

G = sin 2 9 = const. 

associée à la variable cyclique (p. 1 Les solutions (p(t) = po sont des solutions parti- 
culières : mouvements d’un pendule simple dans un plan méridien. Les solutions du 
système précédent ne s’expriment pas en terme de fonctions élémentaires. 

Dans le cas g — 0, les solutions des équations de Lagrange pour un point matériel 
libre sur la sphère sont les géodésiques de la sphère. Il est clair que les méridiens 
(9(t) = 9 0 t + 0o, (p(t) = (fo) sont des géodésiques. Comme l’axe sud-nord a été 

choisi de manière arbitraire, ce sont en fait tous les grands cercles qui constituent 

les géodésiques de la sphère. □ 

1.3.4 Le couplage minimal au champ électromagnétique 

Comment maintenant déduire les équations de Lorentz gouvernant le mouvement 
d’une particule de masse m et de charge électrique q soumise à l’action du champ 
électromagnétique extérieur ? 

Nous savons que ces équations du mouvement ne se déduisent pas de la théorie du 

champ électromagnétique de Maxwell, qu’elles relèvent en quelque sorte d’un autre 

principe : le principe du couplage minimal. Le formalisme lagrangien permet 

une formulation simple et élégante de ce principe de base de l’électrodynamique. 

7. On dit qu’une variable, disons y, est cyclique si dL/dp = 0. Grâce aux équations de 
Lagrange Y impulsion associée 

est automatiquement une intégrale première, p v = const., des équations du mouvement. 
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Soit L 0 = |m||r|| 2 le lagrangien libre de la particule non relativiste. Si l’on 
soumet ce système à l’action du champ électromagnétique (E. B), la prescription du 
couplage minimal consiste à remplacer le lagrangien L 0 par le lagrangien suivant 

(1.3.18) 

qui dépend explicitement du potentiel vecteur A et du potentiel scalaire 0 dont 
dérive le champ électromagnétique selon 


L{r,r,t) = -m\\r\\ 2 + q 


(A(r,t),r) - 0(r ,t) 


E 


— grad0 — 


dA 

~dt 


et B = rotA. (1.3.19) 


Remarque 1.3.7. Nous verrons que les équations de Lagrange pour (1.3.18) ne 
mettent en jeu que le champ électromagnétique et pas les potentiels (non physiques) 
définis à une transformation de jauge près. 


Ecrivons (1.3.18) comme L(r,r,t) = \mY^j= i ( x j ) 2 + q Yl]=i AjX j — (j où les ad 
désignent les composantes (cartésiennes) de la position r G R 3 et ad celles de la 
vitesse r G M 3 ainsi que Aj = dL celles du potentiel (co)vecteur A G C'°°(M 3 x M, M 3 ) 
choisi. 


La i-ème composante de l’impulsion p = dL/di est ainsi 


dL 

Pi = 7777 = mxi + qAi 
ox l 


pour tout i = 1,2,3. En notant d \ := d/dx l l’opérateur de dérivée partielle selon x l 
et d t . := d/dt celui de dérivée partielle par rapport au temps t, on obtient aisément 

^ = q[J2 diA i ij ~ ^ 

3 = 1 


pour tout i = 1,2, 3. 
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Les équations de Lagrange s’écrivent alors 


dpi d L 

dt dx i 


mxi + qÀi - q ^ d % A 3 x j — 8$ 

. 3 = 1 

3 i r 3 

rriXi + q djAi H + d t Ai — q djAjX ■ î — d \<f> 

- 3 = 1 J L 1=1 

3 

+ ç -E Fij x J + d t Ai + di<t> 

. t=i 


où l’on a défini le tenseur magnétique antisymétrique 


F i: j = diAj - djAi 

pour tous i,j = 1,2,3 — on vérifie que F 12 = (rotA) 3 = B 3 , etc., grâce à (1.3.19). 
En d’autres termes la matrice 3x3 antisymétrique (Fj) = ( F VJ ) est donnée par 

F = -j( B) 

où j(B) G L(M 3 ) est l’opérateur linéaire défini par 

j(B)v:=Bxv VvGl 3 

où x désigne le produit vectoriel de M 3 

On obtient enfin les équations de Lagrange mr + q [B x r + d t A + grad0] = 0, 
c’est-à-dire, compte tenu de (1.3.19), 

(1.3.20) 

L’équation différentielle (1.3.20) du second ordre est l’équation de Lorentz dont les 
solutions r (t) constituent les mouvements (on trajectoires) d’une particule chargée 
dans un champ électromagnétique extérieur. 





Chapitre 2 


Les équations de Hamilton 


2.1 Equations canoniques 


Le formalisme lagrangien met enjeu, nous l’avons vu, l’espace “tangent” à l’espa- 
ce de configuration d’un système mécanique, c’est-à-dire l’espace des configurations 
et des vitesses — vecteurs tangents — possibles pour ledit système. Nous avons, 
d’autre part, mis en évidence la notion d’impulsion généralisée qui intervient dans 
la formulation des équations de Lagrange. Dans certains contextes, il peut s’avérer 
judicieux de travailler plutôt dans l’espace dit “cotangent” à notre espace de confi- 
guration, i.e. l’espace des configurations et des impulsions — (co)vecteurs tangents. 
C’est précisément ce que se propose de faire le formalisme hamiltonien dû, entre 
autres, à W.R. Hamilton. 

Soit L( q, q) le lagrangien d’un système (éventuellement soumis à des contraintes 
holonomes) dont l’espace de configuration est cartographié par un système de coor- 
données q = (q 1 , . . . ,q n ) G M C 1" et l’espace tangent au point q décrit par les 
vitesses q = (q 1 , . . . , q n ) G M n . 


Définition 2.1.1. On appelle impulsion généralisée pour le lagrangien L(q, q) 7 
dans le système de coordonnées considéré, le covecteur 



au point q G M , ou encore Pi = dL/dq 1 pour tout i — 1, . . . , n. 


( 2 . 1 . 1 ) 
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Nous voyons que p est alors une fonction de q et q définie par (2.1.1) de sorte 
que l’application 

F l ■ (q,q) ^ (q,p) 

constitue une application (locale) du tangent au cotangent de l’espace de configura- 
tion. Nous supposerons cette application inversible et d’inverse différentiable, pour 
que q puisse aussi être vue comme fonction différentiable de p et de q. 

Nous prouvons maintenant le théorème fondamental suivant : 

Théorème 2.1.2. Le système des équations de Lagrange (1.3.1) est équivalent au 
système de 2 n équations différentielles du premier ordre, les équations de Hamilton , 1 

( 2 . 1 . 2 ) 

pour i — 1, . . . , n où 

(2.1.3) 

est le hamiltonien du système, fonction différentiable H(pi, . . . ,p n , q 1 , . . . , q n ) de 

l’espace des phases . 2 




Démonstration. Calculons la différentielle de L{ q, q) : on obtient alors facilement 
d.L = YÉÏ=i ( dL/dq l )dq l + (dL / dcf)dcf ou encore 


, r dL dL 

d.L = — d.q T ~^~r d.q 

O q O q 


= p • d.q + p • d.q 


en utilisant les équations de Lagrange dL/dq = p (cf. (1.3.1)) et (2.1.1). Une 
intégration par parties donne alors d.L = p ■ dq + d( p • q) — dp ■ q ou encore 

1. On appelle aussi ces équation canoniques. 

2. L’hamiltonien H( p, q) est la transformée de Legendre (2.1.3) du lagrangien L(q, q). L’hamil- 
tonien est défini sur l’espace des phases T*M = (K™)* x M, c’est-à-dire l’espace “cotangent” de 
l’espace de configuration M C K”, alors que le lagrangien est, lui, défini sur le l’espace “tangent” 
TM = K" x M de l’espace de configuration M. 





2.1. EQUATIONS CANONIQUES 


19 


d(L — p • q) = p • dq — dp ■ q que l’on réécrit en définissant, cf. (2.1.3), 

H = p ■ q — L 


comme 

dH = dp ■ q — p ■ dq (2-1-4) 

en mettant en évidence qne l’hamiltonien est une fonction id(q, p). On trouve alors 
aisément q = dH/dp et p = —dH/dq. □ 

Exercice 2.1.3. Calculer le hamiltonien de l’oscillateur harmonique dont le lagran- 
gien est donné dans l’Exercice 1.1.3. En déduire les équations de Hamilton. 

Exercice 2.1.4. Supposons que l’énergie cinétique d’un système mécanique soit une 
fonction T( q, q) et l ’ énergie potentielle donnée par une fonction V (q) . Le lagrangien 
du système est, on le sait, 

L = T — V. 

Trouver la condition à imposer à l’énergie cinétique T pour que le hamiltonien soit 
de la forme suivante 

(2.1.5) 

Démonstration. Indication : prouver, en utilisant la formule d’Euler pour les fonc- 
tions homogènes de degré k, 3 qne T doit être une fonction homogène en la variable q 
de degré k déterminé. □ 


H = T + V 


Remarque 2.1.5. Le Théorème 2.1.2 reste inchangé dans le cas non stationnaire 
où le lagrangien L(q, q, t) dépend explicitement du temps. On trouve dans ce cas que 
les équations de Newton (2.1.2) se complètent par 


dH _ dL 
~dt ~ ~~dt' 


( 2 . 1 . 6 ) 


3. Rappelons qu’une fonction /( r) est homogène de degré k si /( Ar) = A A '/( r ) pour tout A > 0. 
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2 . 1.1 Exercices 

Voici maintenant quelques exercices simples pour illustrer l’usage des équations 
de Hamilton. 

Exercice 2.1.6. Donner le hamiltonien H(pg,9) du pendule simple de longueur I 
dans le champ de pesanteur d'accélération g = const. > 0. Ecrire les équations de 
Hamilton. Conclusion ? 


Exercice 2.1.7. On considère le hamiltonien 


H{ P, r ) = 


(2.1.7) 


fonction de (p, r) G (M 3 \{0}) xl 3 , oùn(r) > 0 une fonction différentiable. Montrer 
que les équations de Hamilton correspondent aux équations de Fermât (1.2.2) pour 
un indice de réfraction variable n( r). 


Démonstration. Indication : introduire la vitesse unitaire u = p/||p|| = dr/ds. □ 


Exercice 2.1.8. Etudions la dynamique hamiltonienne dans un référentiel tournant 
(par exemple lié à un manège) à une vitesse angulaire constante u ) = ue z autour de 
la direction e z . On suppose que ces deux référentiels coïncident au temps t — 0. 

1. Montrer que si q = (x, y, z ) désignent les coordonnées dans le référentiel “fixe” 
et = Q = (. X , Y, Z) celles associées au référentiel tournant, à l'instant t, on a 


x = X cosc ot — Y sin eut (2.1.8) 

y = X sin ut + Y cos ut (2.1.9) 

z = Z (2.1.10) 


2. Exprimer l’énergie cinétique T{ Q,Q) d 'un point matériel libre de masse m. 

3. Si V (Q) désigne le potentiel dans lequel est plongé le point matériel, donner 
l’expression du lagrangien L( Q. Q) du système. 

4- En déduire le hamiltonien H( P.Q). 
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5. Ecrire les équations de Hamilton du système dans le référentiel tournant. ( On 
pourra poser W — X + ïY G C.) 

Démonstration. Réponses : (i) Le lagrangien est donné par 

L{ Q, Q) = X 2 + Ÿ 2 + Z 2 + co 2 (X 2 + Y 2 ) + mw{XŸ - YX) - V(X, Y , Z) 

et l’impulsion canonique donnée par P = dL/dQ = m(X — wY, Y + cjX, Z), (ii) Le 
hamiltonien prend alors la forme 

m p , Q) = A [pt + p; + pt] - ^(xp, - y po + v(x, y, zy 

Enfin (iii) les équations de Hamilton conduisent au système d’équations différentielles 
couplées 


X = 

2uŸ + YX--^~ 
m oX 

(2.1.11) 

Ÿ = 

-2 uX + YY--^ 
m Oï 

(2.1.12) 

Z = 

1 dV 
m d Z 

(2.1.13) 


Remarquons que les deux premières équations différentielles s’écrivent, en introdui- 
sant le changement de coordonnées W = X + ïY suggéré plus haut, 

dV 

mW = —2imujW + muj 2 W — 2^=. (2.1.14) 

dW 

Le premier terme dans le membre de droite de (2.1.14) est la force de Coriolis, le 
second la force centrifuge et le dernier la force extérieure dérivant du potentiel V. □ 

2.1.2 Le couplage minimal au champ électromagnétique 

Donnons ici la version hamiltonienne du principe de couplage minimal à un 
champ électromagnétique, formulé dans le cadre lagrangien dans le Chapitre 1.3.4. 

Nous nous proposons de calculer maintenant l’hamiltonien H( p, r, t) déduit du 
lagrangien (1.3.18) — via la transformation de Legendre (2.1.3) — d’une particule 
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de masse m et de charge électrique q plongée dans une champ électromagnétique 
extérieur (E.B) dérivant, cf. (1.3.19), d’un potentiel vecteur A et d’un potentiel 
scalaire 0. 

Calculons donc le moment conjugué p = dL/dr. Il vient aisément 

p = mv + qA 

de sorte que r = (p — qA) /m. On a alors 
H = p.r — L 

= p • (p — qA )— — L 
m 

1 f 1 1 2 r 1 i 

= P- (p-çA)- - I -m — (p — qA) + q [(A, -(p - qA)) - 0 

111 2 

= (p - ?A) • (p - qA) -m — {p - qA) + q(j) 

m 2 m 

= \\p - qM 2 + qtp- 

2 m 

Si Hq(jp, r) = || p || 2 / (2m) désigne l’hamiltonien libre d’une particule (non relativiste) 
de masse m, l’hamiltonien décrivant les mouvements de cette particule dans un 
champ électrique dérivant d’un potentiel 0 sera donné, cf. (2.1.5), par l’hamiltonien 
standard H(p,r,t) = H 0 ( p,r) + V(r, t) on V = q<j>. Nous venons de prouver que 
le couplage minimal à un champ électromagnétique extérieur s’opère maintenant 
par la prescription pins subtile qui consiste à remplacer H 0 ( p, r) par 

(2.1.15) 

Exercice 2.1.9. Ecrire les équations de Hamilton pour le hamiltonien (2.1.15) et 
retrouver les équations de Lorentz (1.3.20). 

2.2 Crochets de Poisson et transformations cano- 
niques 

Introduisons dans ce chapitre des notions nouvelles et importantes en mécanique 
analytique, à savoir celles de crochet de Poisson et de structure symplectique. Ces 
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notions ont conduit à des généralisations multiples dans le cadre géométrique et 
algébrique relevant des mathématiques et de la physique mathématique. 

2.2.1 Crochets de Poisson 

Considérons une fonction différentiable F : T*M — » M de l’espace des phases 

à valeurs réelles, ce que l’on appelle observable classique F( p, q) en physique. 4 

Comment évolue cette fonction au cours du temps? compte tenu de l’évolution 

temporelle propre du système hamiltonien. Nous avons 

dF d F dp dF dq 

dt dp dt dq dt 

dF dH dF dH 
dp dq + dq <9p 

grâce aux équations de Hamilton (2.1.2). D’où le résultat suivant 


Définition-Théorème 2.2.1. Nous appellerons crochet de Poisson de deux ob- 
servables F et G de l’espace des phases la nouvelle observable 5 




dF dG 
dq 1 dpi 


( 2 . 2 . 1 ) 


de sorte que l’évolution temporelle de toute observable classique F soit gouvernée 
par l’équation différentielle 

( 2 . 2 . 2 ) 

Bien entendu, les équations de Hamilton (2.1.2) sont retrouvées via (2.2.2) car 
Pi = {H, pQ = — dH/dq l et <f = {H, q *} = dH/dpi pour tout i = 1, ... ,n. 

Corollaire 2.2.1. L’hamiltonien H est une constante du mouvement, 

dH 
dt 

4. Nous avons déjà rencontré de telles fonctions, par exemple l’hamiltonien H (énergie), une 
composante pi de l’impulsion d’une particule, une composante q 1 de sa position à un instant donné 
avec (i, j = 1,2,3), etc. 

5. Tout système de coordonnées (pi, . . . ,p n , q 1 , . . . , q n ) dans lequel le crochet de Poisson est de 
la forme (2.2.1) est dit système de coordonnées canoniques. 


(2.2.3) 
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Démonstration. Le crochet de Poisson (2.2.1) est antisymétrique en ses arguments, 
{F, G} = -{G, F}, et donc d.H/dt = {H, H} = 0. □ 

Remarquons que le crochet de Poisson (2.2.1) peut alors s’écrire comme 

{F,G} = ô f G (2.2.4) 


si l’on introduit la dérivation suivante 


F dPi dÉ dÉ dpi 


(2.2.5) 


aussi appelée champ hamiltonien associé à la fonction F de l’espace des phases. 
Le fait que ô F : G h- > {F, G} soit une dérivation implique la propriété importante 
suivante {F, GH} = ô f (FG) = ( Ô F G)H + G5 F H = {F, G}H+G{F, H}. Nous avons 
enfin le 


Théorème 2.2.2. Le crochet de Poisson (2.2.1) est une application bilinéaire 


(F, G) ^ {F, G} 


jouissant des propriétés suivantes 

1. {F,G} = — {G, F} (antisymétrie), 

2. {F, GH} = {F, G}H + G{F, H} (règle de Leibniz), 

3. {F, {G, H}} + {G, [H, F}} + {H, {F, G}} = 0 (identité de Jacobi). 

Exercice 2.2.3. Prouver l’identité de Jacobi (on se restreindra au cas n — 1). 

Exercice 2.2.4. Soient (p, q) les coordonnées canoniques de T*M n . Calculer les 

crochets de Poisson mutuels des composantes pi et q J , pour i, j = 1, . . . ,n. Montrer 

que l’espace vectoriel réel f) n de dimension 2 n + 1 engendré par pi, . . . ,p n , q 1 , ■ ■ ■ , q n 

et la fonction constante 1 est stable sous le crochet de Poisson. 6 

6. On appelle un tel espace algèbre de Lie ; ici f) n est une algèbre de Lie associée au nom de 
Werner Heisenberg, physicien allemand (1901-1976). 
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Exercice 2.2.5. Considérons l’espace des phases T*R 3 et désignons par L = pxqle 
moment angulaire, (i) Calculer pour i,j = 1,2,3 les crochets de Poisson { L i} Lj } des 
composantes du moment angulaire, (ii) Un hamiltonien H( p, q) est invariant sous 
rotations euclidiennes ssi { Li,H } = 0 quel que soit i — 1,2,3. A quelle condition 
l’hamiltonien H = ||p|| 2 /(2m) + E(q) d’une particule de masse m plongée dans un 
potentiel V (q) est-il invariant sous rotations ? 


2.2.2 Structure symplectique 


Nous introduisons dans cette Section une notion nouvelle, celle de forme symplec- 
tique canonique de l’espace des phases T*R n ; cette notion, qui relève de la géométrie 
différentielle à un niveau plus avancé, sera néanmoins traitée brièvement afin d’in- 
troduire à celle de transformations canoniques : changements de coordonnées de 
l’espace des phases qui laissent les équations de Hamilton invariantes, et qui donc 
invarient les lois de la mécanique. 

Réécrivons l’expression (2.1.4) de manière tout à fait équivalente comme 

^ • 6q - ôp • ^ = -ÔH (2.2.6) 

dt M v dt v ' 


pour tous vecteurs (dp, dq) de l’espace des phases; en effet, SH = (dH/d p) • dp + 
(dH/dq) ■ dq représente la dérivée de l’hamiltonien H dans la direction (dp,dq). 
L’expression précédente restitue bien les équations de Hamilton (2.1.2), à savoir 

( dp dH 


dt 

dq 

dt 


dq 

dH 


(2.2.7) 


Définition-Théorème 2.2.2. Définissons, en tout point x = (p, q) de l’espace des 
phases T*M. n = (M n )* x M n l’expression suivante 


u)(Sx, ô'x ) = dp • d'q — <Lp ■ <5q. 


( 2 . 2 . 8 ) 


(i) Cette expression dépend bilinéairement des vecteurs Sx, ô'x G T*R n et est anti 
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symétrique : u u(ôx,ô'x) = —lu(ô'x, ôx). (ii) On a de plus 

ou{ôx , S'x) = ôa(ô'x) — ô'a(ôx) où a(ôx) = p • ùq (2.2.9) 

et (iii) kercu = {0}. Une telle application lu : T*R n x T*R n — > R. est appelée forme 
symplectique canonique de l ’espace des phases. 7 

Démonstration. On a clairement l’identité lu(ôx,ô'x) + lu(ô'x, ôx) = 0 et de plus 
Lu(ôx,aô'x + bô"x) = acu(ôx, ô'x) + bu(ôx,S ,, x ) pour tous a, b e M; d’où le (i). 
L’expression (2.2.9) reproduit simplement la définition (2.2.8), ce qui justifie (ii). 
Enfin, Sx G keru; ssi lu(ôx,ô'x) = 0 pour tout ô'x ; mais, cette dernière condition 
implique à la fois ùp = 0 et ùq = 0, c’est-à-dire ôx = 0, d’où (iii). □ 

Soit H(x) un hamiltonien, le champ de vecteurs x h- > Ôhx défini par (2.2.6), 
c’est-à-dire 

u(ô H x,ô'x) = -Ô'H (2.2.10) 

est le champ hamiltonien associé à H — déjà introduit en (2.2.5) ; cette dernière 
définition est maintenant intrinsèque (car elle ne met en jeu que la fonction H et la 
forme symplectique lu mais pas un système de coordonnées particulier). 

Proposition 2.2.6. Les équations de Hamilton associées à un hamiltonien H(x) 
prennent la forme d’une équation différentielle ordinaire du premier ordre, à savoir : 

f-=S H x ( 2 , 2 , 11 ) 

dt 

avec la définition (2.2.10) du champ hamiltonien Ôh ; le crochet de Poisson de deux 
observables F et G quelconques retient alors la forme suivante : { F , G} = ÔfG. 

7. Il est possible, mais non obligatoire, de reformuler les résultats précédents en termes de formes 
différentielles : lu = dpi A dq 1 est une 2-forme inversible et fermée de T* K" ; elle est, en fait 
exacte puisque u = da avec a = Y)'i=i PNC ■ 
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2.2.3 Transformations canoniques 

Définition intrinsèque 


Les équations de Hamilton (2.2.7) on été écrites dans un système de coordon- 
nées (p, q) particulier de l’espace des phases T*R n . Est-il possible de déterminer 
un changement de coordonnées / : T*R n — > T*W n tel que (P. Q) = /(p, q) et tel 
que les équations de Hamilton retiennent la même forme dans les deux systèmes de 
coordonnées ? à savoir 


d h 

P = -àq 
d h 

q = + âp 


et 


P = 


Q — d 


dH 

dQ 

dH 

dP 


( 2 . 2 . 12 ) 


avec l’expression suivante de l’hamiltonien dans les deux systèmes de coordonnées 


h(p, q) = H(P. Q). (2.2.13) 

Lin tel changement de coordonnées x i— *• X — f{x) doit être différentiable et (locale- 
ment) inversible, x = / _1 (X), et préserver la forme symplectique (2.2.8), c’est-à-dire 

u(ôx, S'x) = u(ÔX, S’X) (2.2.14) 

pour préserver la forme des équations de Hamilton associées au hamiltonien (2.2.13). 

Définition 2.2.7. On appelle transformation canonique 8 toute application dif- 
férentiable /:ih X de T*R n qui invarie la forme symplectique u selon (2.2.14). 

Exercice 2.2.8. A quelle condition une transformation linéaire A : x i— >• X — Ax 
de M 2 est-elle une transformation canonique ? 

Exercice 2.2.9. Soit q i— >• Q(q) une fonction monotone, Q'(q ) 0, d’une variable 

réelle. Vérifier que l’application f : M 2 — ■> M 2 : ( p , q) > (P = p/Q'{q ), Q{q)) est une 
transformation canonique. 9 


8. On dit aussi symplectomorphisme. 

9. Les transformations canoniques introduites dans cet exercice se généralisent au cas d’un 
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Démonstration. On a P8Q = ( p/Q'(q))SQ(q ) = {p/ Q' (q))Q' {q)8q = p8q. Donc 
8P8'Q - 8'P8Q = SpS'q - 5'pôq , d’où (2.2.14). □ 


Introduisons maintenant la notion utile de fonction génératrice d’une transfor- 
mation canonique. 

Lemme 2.2.10. Soit x i— >• A" une transformation de T*R n telle que 


a(8X) — a(ôx) + SF (2.2.15) 

pour une certaine fonction différentiable F : T*R n — > M de l’espace des phases . 10 
Cette transformation est automatiquement une transformation canonique; on ap- 
pelle F fonction génératrice de cette transformation canonique. 

Démonstration. On a uj{8X, 8'X) = 8a(8'X ) — 8'a{8X) = <5(a(<50é)) — ù'(a(ÙA^)) — 
a(66'X - 8' SX). Alors u(SX, 8'X) = u(8x, 8'x ) + {88' F - 8' 8F) - {88' F - 8' 8F) = 
u{8x,8'x). □ 

Nous pouvons donc écrire, en particulier, 

P -8Q-p-8q = 8F (2.2.16) 


ce qui implique que la fonction génératrice est, ici, une fonction F{ Q, q) telle que 


dF 

dQ 


& 



(2.2.17) 


Exercice 2.2.11. Trouver la transformation canonique définie par la fonction génératrice 


suivante :F{ Q,q) = (Q,q). 


Démonstration. On trouve, grâce à (2.2.17), P = q et p = — Q. Il vient donc 

(P,Q) = (q, -p) qui est bien une transformation canonique. 11 □ 

nombre arbitraire n de degrés de liberté. En effet, soit Q : q i— > Q = Q(q) : K” — » K™ une 
transformation différentiable et d’inverse différentiable, on montre (exercice !) que la transformation 
de T*R n suivante / : (p, q) (P, Q) où 



est une transformation canonique. 

10. On supposera F au moins deux fois différentiable. 

11. En effet, on a <5P • Ô'Q — <5'P • <5Q = <îq • <7(~P) — ù'q • ô{— p) = ép • ù'q — é'p • ùq. 
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Exemple de l’oscillateur harmonique 

Traitons ici un cas exemplaire, celui de l’oscillateur harmonique à un degré de 
liberté caractérisé par une masse m > 0 et une pulsation propre u > 0. 


L’hamiltonien du système est, on le sait, 

u n P 2 , mu 2 2 

h[p ' q)= 2Ïi + — q 

et la fonction génératrice envisagée de la forme suivante 

F(Q, q) = — — g cotgQ. 
On trouve aisément, en utilisant (2.2.17), 


mu 2 


2 sin 2 Q 


q 2 > 0 & p = muqcotgQ, 


c’est-à-dire 


p = V 2rnPu cos Q & q — \ sin Q, 

V mu 


formule qui définit la transformation canonique inverse : ( P,Q ) i— >• ( p,q ). 
L’hamiltonien prend alors la forme suivante 

r i i 2 

^ r T 2 TïlÜJ ^ / 2.P 

h(p,q) = V 2mPu cos Q H \/ sinQ 

2m L J 2 V mcn 


= u P cos 2 Q + cnP sin 2 Q 


= uP. 


D’où, en définitive, cf. (2.2.18), 


(2.2.18) 


(2.2.19) 


( 2 . 2 . 20 ) 


h(p, q) — H (P, Q) — uP . 


( 2 . 2 . 21 ) 


La forme symplectique étant invariante par construction, les équations de Hamilton 
dans les nouvelles coordonnées donnent alors une forme très simple aux équations 


du mouvement 




( 2 . 2 . 22 ) 
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qui permet une intégration immédiate 

p(t)=P( 0) & Q(t) = ut + Q(0). (2.2.23) 

Remarquons, pour terminer, que P = H/ui est une constante homogène à une 
énergie/fréquence et représente l’action du système tandis que Q représente la 
phase, ou angle Q(t) = 9(t), du système puisque, grâce à (2.2.20) et (2.2.23), 

q(t) = sin M + W) 

V mu: z 

représente bien la trajectoire de l’oscillateur harmonique déterminée par les valeurs 
initiales des variables action-angle. 



Chapitre 3 


Mécanique des systèmes en 
repères mobiles 

3.1 Le groupe euclidien 

3.1.1 Espace euclidien 

On appelle espace affine associé à un espace vectoriel (réel) V n de dimension 
n un ensemble A n muni d’une action (i) libre et (ii) transitive du groupe additif 
( V n , +). Ceci signifie que l’on s’est donné une loi A n x V n — > A n : (M, v) i — >• M + v 
vérifiant (M + v) + w = M + (v + w) et telle que 

- (i) M + v = M y» v = 0, 

- (ii) si M, N G A n il existe un unique vecteur v G V n tel que N = M + v. 

On note 

MN = N ~ M = v (3.1.1) 

le vecteur d’origine M et d’extrémité N associé au bipoint (. M,N ). 

Nous voyons donc que l’on peut identifier A n et V n dès qu’un point O G i" a 
été choisi : l’isomorphisme (affine !) A n — > V n est alors donné par 

M ^ OM. 

La définition suivante nous apprend à repérer les points d’un espace affine. 
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Définition 3.1.1. On appelle repère (affine) d’un espace affine (A n , V n ) un couple 
IZ = (O, S) formé d’une “origine” O G A n et d’une base S = (ei . . . e„) de V n A 
Les coordonnées r = (a: 1 , . . . ,x n ) G M n d’un point M , dans ce repère, sont définies 
par 

M = IZr 
= O + Sr 

n 

= o + Gj x i . 

i = 1 

Exercice 3.1.2. Montrer que les changements de repères affines forment un groupe, 
le groupe affine composé des matrices (n + 1) x (n + 1) de la forme 

°=( 1 ?) < 31 - 2) 
où A G GL(n, t), le groupe multiplicatif des matrices réelles n x n inversibles, 
et b G M n . 

Définition 3.1.3. On appelle produit scalaire euclidien d’un espace vectoriel V n 
toute application bilinéaire symétrique non dégénérée positive : 

g : V n x V n -> M 

c’est-à-dire vérifiant 

1. g(X ivi + À 2 v 2 , w) = Ai gfvi, w) + A 2 g(v 2 , w) 

£ g(v, w) = g( w,v) 

5. ^(v, w) = 0 pour tout w G K" ssî v = 0 
|| v|| 2 := g(v,v) > 0 [et ||v|| = 0 ssi v = 0/. 
pour tous Ai, A 2 G M et vi, v 2 , v, w G V n . 

Définition 3.1.4. Nous dirons qu’une base S = (e i . . . e n ) est orthonormée si 

g(e i ,e j ) = 8 ij (3.1.3) 

pour tous i,j = 1, . . . , n. 2 


1. Une base est un isomorphisme linéaire S : K” — > V n . 

2. On désigne par 5 lJ le symbole de Kronecker, égal à 1 si * = j et à 0 sinon. 
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On note g(v,w) = (v,w) le produit scalaire euclidien canonique de M" donné 

par 

n 

(v,w)=vw = ^v i w i (3.1.4) 

i= 1 

en désignant par une barre la transposition. 3 

Nous pouvons maintenant donner la définition générale d’un espace euclidien. 

Définition 3.1.5. On appelle espace euclidien tout espace affine E n dont V espace 
vectoriel associé est un espace euclidien ( V n ,g ). 

Si une base S orthonormée de V n est donnée on notera 

g(SM,ô'M) = (6r,ô'r) (3.1.5) 

le produit scalaire de deux vecteurs JM, J'M G V n d’origine M G E n et de compo- 
santes ôr, S ' r G M n dans cette base. 

3.1.2 Isométries euclidiennes 

Le groupe orthogonal 

Considérons l’ensemble des matrices n x n réelles, A, inversibles, dont l’inverse 
est égal à la transposée, c’est-à dire telles que 

A” 1 = 2. (3.1.6) 

Ces matrices forme un groupe. En effet, A = 1 (la matrice identité) vérifie bien (3.1.6) 

et sera clairement l’élément neutre. De plus, si A et B sont deux telles matrices, il 

3. La transposée A d’une matrice A h m lignes et n colonnes est la matrice à n lignes et m 
colonnes obtenue en échangeant lignes et colonnes de A , i.e. (Af = A j pour tous i = 1, ... ,n et 
j = l,...,m. 

En particulier le transposé d’un vecteur 



est le covecteur 


v= {v 1 ...v n ) G (K”)*. 
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en va de même de leur produit (matriciel) : (AB)- 1 = R^A -1 = B A = AB. Enfin 
l’inverse B = A -1 de A vérifie bien B - 1 = A = A -1 = B. Remarquons que les 
matrices (3.1.6) sont telles que det(A) 2 = 1. 

Définition 3.1.6. On appelle groupe orthogonal le groupe 

O (n) = {A G L(R n ) | AA = 1} 

Le sous-groupe 

SO (n) = {A G O (n) | det(A) = 1} (3.1.7) 

est le groupe des rotations euclidiennes IRA (ou groupe Spécial Orthogonal). 

Les transformations de IRA données par les matrices orthogonales préservent le 
produit scalaire euclidien ; on a bien (Au, Av) = ( AAu, v) = (u, v) pour tous 
u, v G IRA. Le groupe orthogonal est un groupe d'isométries. 4 

Rappelons que l’espace IRA est orienté par le choix d’une forme volume ; un choix 
traditionnel est donné, pour n vecteurs v 1; . . . , v n G IRA, par le volume suivant 

vol(vi, . . . , v n ) = det(vi . . . v n ) (3.1.8) 

c’est-à-dire par le déterminant de la matrice n x n dont les colonnes sont constituées 
des n vecteurs considérés. 

Remarque 3.1.7. Rappelons que dans le cas n — 3 on a la relation suivante 

vol(vi, v 2 , v 3 ) = (vi, v 2 x v 3 ) (3.1.9) 

entre forme volume, produit scalaire et produit vectoriel. 

Proposition 3.1.8. Toute matrice de rotation A G SO(3) est de la forme 

A = (u v w) (3.1.10) 

cm u, v G M 3 sont des vecteurs unitaires, ||u|| = ||v|| = 1, et orthogonaux, (u, v) = 0, 
et où w = u x v. 


4. Transformations qui préservent la métrique (alias le produit scalaire). 
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Le cas n — 3 est, par ailleurs, d’une importance toute particulière en mécanique 
du solide ; illustrons-le par l’exercice suivant. 

Exercice 3.1.9. Désigner, parmi les matrices suivantes, celles qui sont des rotations 
euclidiennes et celles qui sont des matrices orthogonales : 

( cos 9 — sin 9 0 \ / cos 9 — sin 9 0 \ 

sin 9 cos 9 0 J , D — J sin 9 cos 9 0 J , 

0 01/ \ 0 0-1/ 

/ -1 0 0 \ / 1 0 0 \ 

P = 0-1 0 , S= 0 1 0 . 

\ 0 0 - 1 / \ 0 0 -1 / 

Les rotations, matrices orthogonales A de déterminant positif, préservent 
l’orientation de l’espace : vol(Au, Av, Aw) = det(A)vol(u, v, w) = vol(u, v,w) 
puisque det(A) = 1. 

Les symétries, matrices orthogonales S de déterminant négatif, renversent 
l’orientation : vol(Su, SV, Sw) = — vol(u, v,w) puisque det(S') = —1. 

Le groupe des translations 

Si l’on considère un espace euclidien E n modelé sur M n et muni du produit 
scalaire g, alors les translations (3.1.1), 

M i-> M* = M + b 

avec b G M n , préservent le produit scalaire (3.1.5) car <5M* = <5M puisque b est un 
vecteur constant, c’est-à-dire g(àM, PM) = ^(àM*, PM*) ou encore, si àr (resp. àr*) 
représente <5M (resp. PM*) dans une base orthonormée, 

(Sr, Pr) = (<5r*, Pr*). (3.1.11) 


Les isométries euclidiennes 

Nous pouvons conclure que les transformations d’un espace euclidien ( E n ,g ) de 
la forme 


r i— ► r* = Ar + b 


ou 


A G O (n) & bGP 


(3.1.12) 
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sont des isométries, i.e. préservent le produit scalaire g. Ce sont des transformations 
affines (3.1.2) particulières. Nous admettrons le résultat suivant : 


Théorème 3.1.10. Les isométries d’un espace euclidien ( E n ,g ) sont constituées 
des transformations (3.1.12) qui forment un groupe appelé le groupe euclidien : 


E(n) 


A b 

0 1 


A e 


O (n), b G M. 



(3.1.13) 


Le sous-groupe des transformations euclidiennes qui préservent l’orientation est ap- 
pelé groupe spécial euclidien : 


SE (n) 


A b 

0 1 


A g 


SO (n), b G R. 



(3.1.14) 


Exercice 3.1.11. Expliciter la loi de groupe du groupe euclidien : (0) donner l’élément 
neutre, (i) calculer le composé [A!', b") de {A, b) et (A',b') et (ii) trouver {A, b) -1 . 


3.2 Changements de référentiels non inertiels 

3.2.1 Prolégomènes 

L’espace-temps non relativiste est, on le sait, un espace galiléen, i.e. un espace 
affine, A 4 , muni (i) d’une fonction temps absolu t : A 4 — > A 1 et (ii) d’une structure 
d’espace euclidien E ?> = ( A 3 ,g t ) sur chaque espace instantané t = const. 

Dans un référentiel galiléen IZ où un événement d’espace-temps est représenté 
par sa position et sa date, 

m = k (0 

avec r G l 3 et! G I, on a 

<7t(<5M, ù'M) = (ôr, ô'r) si St — S't — 0 (3.2.1) 

et dt désigne la différentielle donnant l’intervalle de temps entre deux événements. 

Les symétries galiléennes, c’est-à-dire les transformations différentiables d’espace- 
temps M i— > M* préservant la structure “métrique” galiléenne (g t , dt) et l’orientation 
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spatiale des espaces t = const. sont donc données par 


r* = A(t)r + b (t) 
t* = t + e 


(3.2.2) 


où, cf. (3.1.13), les matrices A(t) G S0(3) et les vecteurs b(f) G M 3 peuvent mainte- 
nant, grâce à (3.2.1), dépendre arbitrairement du temps ; d’autre part e G M désigne 
une translation temporelle (solution générale de l’équation dt* = dt). 


Définition-Théorème 3.2.1. Les transformations (3.2.2) forment un groupe (de 
dimension infinie) ; le sous-groupe formé des transformations préservant les espaces 
instantanés t = const., i.e. 

+ (3.2.3) 

est donc le groupe des fonctions ( A , b) G C 00 ^, SE(3)) appelé groupe de Coriolis 
ou groupe des transformations non inertielles. 



Enonçons maintenant un résultat fort utile. 


Lemme 3.2.1. Soit A(t) G SO(n) une matrice de rotation dépendant différentiable- 
ment du temps t. La matrice Z(t) = À^Alf)^ 1 est alors antisymétrique, c’est-à-dire 
Z(t) + Z (t) = 0. Dans le cas n = 3, on a 


À^Ait)- 1 = j{u(t)) 


où u }(t) G M 3 est appelé vecteur instantané de rotation et 


/ 0 —en 3 en 2 \ 

[ lu 3 0 —en 1 ] 
V -en 2 en 1 0 J 



(3.2.4) 


(3.2.5) 


Démonstration. Puisque AA = AA = 1, on a d(AA)/dt = 0 et on trouve, par 
conséquent, ( dA/dt)A + A(dA/dt ) = (cL4/ dt)A~ l + {dA/ dt)A~ l = Z + Z = 0 . Dans 
le cas n— 3, on vérifie directement que les matrices 3x3 antisymétriques sont bien 
de 1a. forme générale (3.2.5). □ 
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Exercice 3.2.2. Déterminer le vecteur instantané de rotation lo de la matrice A de 
l’Exercice 3.1.9 dépendant du temps via une fonction 9{t). 

Exercice 3.2.3. Vérifier que le produit vectoriel de lo,lo' G M 3 est donné par 

(3.2.6) 

Exercice 3.2.4. Montrer que le double produit vectoriel 

lo x (lo' x lo") =lo'(lo,lo") — <o " (lo , lo') (3.2.7) 

de trois vecteurs lo,lo',lo" G M 3 est donné par 

j(io)j(io') = lo'lô-lôlo 1 . (3.2.8) 

Exercice 3.2.5. Vérifier l’identité suivante 

j(Lo)j{Lo') - j(io')j(io) =j(iox LO ') (3.2.9) 

pour tous lo,lo' G M 3 . 

Nous utiliserons enfin le résultat précieux : 

Proposition 3.2.6. Soient A G SO(3) et lo G M 3 , on a alors 

Aj^A- 1 = j(Aw). (3.2.10) 

Démonstration. Nous avons, pour toute matrice A G 0(3) et tous u, v,lo G M 3 , 

(fiu, x fiv) = vol (Au, A*;, Av) 

= det(A)vol(u, lo, v) 

= det(4)(u,w x v) 

et, grâce à (3.1.6) et (3.1.9), on obtient (u, A _1 j(Au;)Av)) = det(A)(u, j(co)v) pour 
tous u, v G M 3 , c’est-à-dire 

A~ 1 j(Alo)A = det(A)j(co) si A G 0(3). (3.2.11) 



Le fait que det(A) = 1 si A G SO(3) achève la preuve. 


□ 
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3.2.2 Considérations mécanistes 

Introduisons maintenant ces transformations sous une forme moins abstraite qui 
utilise un solide de référence “fixe” dans l’espace tridimensionnel M 3 associé à un 
référentiel galiléen ; notons O le point du solide correspondant à l’origine r = 0 et 
(ei, e 2 , e 3 ) la base canonique de M 3 ; le référentiel euclidien IZ = (0, ei, e 2 , e 3 ) sert à 
repérer les points AI de l’espace et on note r = OM le rayon vecteur d’origine O et 
d’extrémité M. Attachons maintenant un référentiel euclidien 7 Z' = ( O e) , e 2 , e 3 ) à 
un autre solide évoluant dans l’espace au cours du temps t. 

Comment repérer le même point M relativement au référentiel mobile 1Z' ? 

Ecrivons la relation de Chasles OM = OO' + O'M, ou encore r = b + r' en 
notant b = O O' le changement d’origine et r' = r re i = O M la position relative 
du point AI par rapport à O'. On désigne traditionnellement par R les coordonnées 
du vecteur r re i dans la base (e 3 e 2 e 3 ) = (ei e 2 e 3 )A déduite de la base originelle 
(canonique) par une rotation A. 5 On a alors 

r re i = AH. (3.2.12) 

Soulignons que rotations et translations sont en réalité des fonctions A{t) G SO(3) 
et b(f) G M 3 du temps t (fonctions que nous supposerons différentiables) puisque le 
deuxième solide de référence est mobile au cours du temps. 

En résumé, les coordonnées, r, du point AI dans le repère “fixe” et ses coor- 
données, R, dans le repère mobile sont reliées par la formule simple 

r = A(t)R + b(t) (3.2.13) 

qui n’est rien d’autre que l’expression de la transformation de Coriolis (3.2.3) — à 
un changement de notation près. 

Appliquons maintenant les résultats précédents au calcul de loi de transformation 
de la vitesse et de l’accélération sous un changement de référentiel correspondant à 
une transformation de Coriolis (3.2.13). 

5. Nous identifierons la matrice de rotation A à la base mobile (e( e' 2 e 3 ) puisque (ej e 2 e 3 ) = 1 
est la base canonique de K 3 . 
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3.2.3 Loi de transformation de la vitesse 

Si un point matériel se déplace dans le repère “fixe” sur une trajectoire Mit ), sa 
vitesse “absolue” est, bien sûr, v(= v abs ) = r. Quid de sa vitesse relative 

v rel = /IR, (3.2.14) 

i.e. de sa vitesse relativement au repère mobile ? 

En différentiant (3.2.13) par rapport au temps, il vient 

r = ÀR + ÆR + b 

= ÀA^ARY AR + b 

= j(u)AR + AR + b (3.2.15) 

= u> x r re i + dR + b 

grâce à la définition (3.2.4) du vecteur instantané de rotation dont on note u:(t) les 
composantes dans le repère fixe. 

Proposition 3.2.7. La vitesse absolue, v, d’un point matériel est reliée à sa vitesse 
relative, v re i (définie par (3.2. If)), par l’expression suivante 

(3.2.16) 

Exercice 3.2.8. Montrer que la loi (3.2.16) de transformation de la vitesse se lit 
de la manière suivante 

v = AfR + flxR + A^bj (3.2.17) 

en terme des composantes dans le repère mobile ; on a posé 

lo = ATI (3.2.18) 

pour définir les composantes Cl du vecteur instantané de rotation dans le repère 
mobile. 
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3.2.4 Loi de transformation de l’accélération 

Nous pouvons maintenant tirer avantage des résultats précédents concernant la 
loi de transformation de la vitesse lors d’un changement de référentiel non inertiel 
pour déterminer celle de l’accélération. 

L’accélération absolue du point matériel est, bien entendu, a(= a a b s ) = r. Comme 
en (3.2.14), définissons naturellement l’accélération relative du point par rapport au 
repère mobile selon 

a re i = AR. (3.2.19) 

Proposition 3.2.9. L’accélération absolue, a, d’un point matériel est reliée à son 
accélération relative, a rel ( définie par (3.2.19)), par l’expression suivante 

(3.2.20) 

Démonstration. En différentiant (3.2.15) par rapport au temps t, il vient 

r = j(cb)4R + j(w)4R + j(w)4R + ÆR + AR + b 

= j(ù})AH + j(u>)j(u>)AR, + j(u>)AR + j(u>)AR, + AR + b 
= lo x 4R oj x (cü x 4R) T x 4R, 4R ~h b 
= c i) x r re i + U) x (u) x r re i) + 2u x v re i + a re i + b 

grâce aux définitions ci-dessus de la position, vitesse et accélération relatives. □ 

Exercice 3.2.10. Montrer que la loi (3.2.20) de transformation de l’accélération se 
lit de la manière suivante 

a = A ^R + 2f2 xR + f2xR + f2x(f2xR) + 4 _1 b) (3.2.21) 

en terme des composantes dans le repère mobile. 



Démonstration. Indication : utiliser (3.2.17). 


□ 
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3.2.5 Forces inertielles : introduction générale 


Considérons un point matériel de masse m, soumis, à l’action d’une force dont 
l’expression est f(r,r,f) dans un référentiel inertiel donné. Comment formuler les 
équations dn mouvement du système dans un référentiel non inertiel déduit du 
précédent par une transformation (3.2.13) ? 11 suffit simplement d’utiliser le résultat 
fondamental donné par l’expression (3.2.20) ! 

La deuxième loi de Newton 

ma = f (3.2.22) 

peut se réécrire, grâce à (3.2.20), comme 

f = m ^a re i + 2u? x v re i + Ce x r re i + w x (w x r re i) + b^ , 

ou encore comme 

ma re i = f — 2mu> x v re i — mû ; x r re i — mu? x (u? x r re i) — mb (3.2.23) 

en mettant en évidence de nouvelles forces intervenant dans la formulation des 
équations du mouvement en référentiel accéléré, les forces d’inertie. 6 

Définition-Théorème 3.2.2. Les équations gouvernant le mouvement d’une parti- 
cule de masse m se formulent comme suit dans un référentiel non inertiel 

ma r el = l"rel = f T fcor T feentrif H - fent (3.2.24) 

où les forces d ’inertie sont respectivement 1 


fCor 

= —2 mu? x v re i 

(force de Coriolis) 


feentrif 

= —mu? x (u? x r re i) 

(force centrifuge) 

(3.2.25) 

f 

A ent 

= —mû? x r re i — mb 

(force d'entrainement) 



6. On rencontre aussi le terme “forces fictives” dans la littérature ancienne; cette terminologie 
est trompeuse car les forces d’inertie sont des forces bien réelles comme chacun peut en faire 
l’expérience dans la vie courante. 

7. Attention ! Les formules (3.2.25) donnent bien les composantes des forces inertielles, mais 
exprimées dans le repère . . . fixe ! 
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Déduisons enfin de (3.2.23) une formulation alternative des équations du mou- 
vement dans un référentiel accéléré arbitraire. 

Théorème 3.2.11. Les équations gouvernant le mouvement d’une particule de 
masse m se formulent comme suit dans un référentiel non inertiel 

(3.2.26) 

en terme des grandeurs cinématiques et dynamiques relatives, où. F défini par 

f(r, r, t) — A(t) F (R, R, t) (3.2.27) 

représente les composantes de la force extérieure dans le référentiel mobile. 

Exercice 3.2.12. Désignons par g le champ de gravitation newtonien de la terre. 
Un pendule est au repos par rapport à la terre, au voisinage du sol; déterminer 
son accélération relative a re i en fonction de g, de la vitesse angulaire u ; de la terre 
lors de son mouvement diurne et de la position r re i = const. du pendule. En déduire 
l’expression g(\) de l’intensité de l’accélération de la pesanteur sur terre en fonction 
de la latitude X, de la vitesse angulaire en = ||en||, du rayon R de la terre et de la 
valeur de go = g(± 7t/2) aux pôles. Af.B. Négliger les termes 0(cn 4 ). 

Démonstration. Puisque v re i = 0 et en = const. (la vitesse angulaire de la terre par 
rapport à un référentiel “fixe” copernicicn est constante) et b = 0, (3.2.23) nous 
donne a rel = g - en x (en x r rel ). Alors g 2 = ||a re i|| 2 = ||g|| 2 - 2(g,en x (en x r rel )) en 
négligeant des termes d’ordre 4 en en = 27r/(24h) = 7,27 10 -5 s -1 . On obtient donc 
g 2 = gl~ 2(g, en) (en, r re i) + 2cn 2 (g,r re i) grâce à (3.2.7) et au fait que go = ||g|| aux 
pôles où r rel //en. Le champ g = — kr re \/R 3 étant central (ici k = const. > 0), on a 
(g, r re i) = —goR au voisinage dn sol. Mais (en, r re i) = uRcosfïï /2 — À) = enRsinÀ 
et (g, en) = g 0 u> cos(7t/2 + À) = — e/o^sinÀ. Alors g 2 = g\ + 2goOJ 2 Rshf 2 X — 2g 0 uj 2 R 
entraîne g 2 = g% — 2g 0 (v 2 Rcos 2 À, c’est-à-dire 

g(X) = g 0 — tu 2 R cos 2 À. (3.2.28) 


mR = F — 2md x R — rnÙ x R — mfi x (12 x R) — rriA 1 b 


On vérifie que g(±n/2) = g 0 aux pôles et g m \ n = ^(0) = go — u 2 R à l’équateur. □ 
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3.2.6 Exemple : chute libre et déviation vers l’est 


Appliquons les résultats généraux précédents au cas de la chute libre dans les 
référentiels accélérés. 

Nous présenterons ici le calcul de la trajectoire d’un point matériel en chute libre 
par rapport à un référentiel terrestre en rotation uniforme, autour de l’axe des pôles, 
de vitesse angulaire u par rapport à un référentiel inertiel copernic-ien fixe. 

Nous effectuerons l’intégration (approchée) des équations du mouvement de ce 
point matériel dans le repère mobile attaché à la terre défini par la matrice 8 


A = (e 0 e p e r ) 


où (r, 9, ip) est un système de coordonnées polaires adapté au problème : le pôle Nord 
correspond à la colatitude 9 = 0 et le pôle Sud à 9 = n. 

Le vecteur instantané de rotation terrestre est donné par u) = u e 3 = ATI avec 
u> = const. > 0, c’est-à-dire 

/ — sin 9 \ 

ft = u;[ 0 . (3.2.29) 

y cos 9 J 


Notons que le champ de gravitation est central et . . . attractif, i.e. 



(3.2.30) 


où g = const. > 0 désigne l’intensité de l’accélération de la pesanteur au voisinage 
du sol. 

Les équations du mouvement (3.2.26) s’écrivent alors 


R = G — 2 fi x R — x (12 x R) 


(3.2.31) 


puisque fl = const. et b = 0. 

8. La matrice de changement de base A G SO(3) est identifiée à la base mobile — orthonormée 
directe — attachée à la terre. 
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Nous nous proposons d’intégrer (3.2.31) en négligeant 9 les termes en 0(cu 2 ). 
Les conditions initiales choisies seront R(0) = Ro et R(0) = 0, i.e. on étudie la 
trajectoire d’un point matériel lâché sans vitesse initiale d’un point R 0 . Les équations 
du mouvement seront donc approximées par 


R = G — 212 x R. 


(3.2.32) 


Posons a priori R(t) = |G t 2 + R 0 + Ri(t) où Ri(t) désigne la petite déviation 
de trajectoire recherchée. Puisque ||Ri|| = O(ca), on trouve que (3.2.32) entraîne 
Ri = 2G x 12 1 modulo des termes du second ordre en u. Une intégration élémentaire 
donne Ri = G x O t 2 puisque fl et G sont des vecteurs constants et Ri(0) = 0. 
Il vient alors Ri(t) = G x flt 3 /3 puisque Ri(0) = 0. On obtient finalement la 
trajectoire suivante 

R(t) = ^G x 12t 3 + ^Gt 2 + R 0 . (3.2.33) 

O Zj 


En introduisant la latitude À = 7t/2 — 6 du point considéré réécrivons (3.2.33) 
comme 


/ 


R(t) - Ro = 


° ^ 
geo cos À — 

t 2 


(3.2.34) 


V -9 2 / 

pour mettre en évidence une déviation vers l’est dans l’hémisphère nord 10 par 
rapport à la verticale lors de la chute libre. 


Exercice 3.2.13. Déterminer la déviation vers l’est d’un point matériel lâché sans 
vitesse initiale du haut de la tour Eiffel, c’est-à-dire d’une altitude H = 275m. On 
prendra À = 49° et g = 10 m s~ 2 . 

Démonstration. Posons R = (X, Y, Z) pour obtenir, grâce à la loi horaire (3.2.34), 
Z — Z 0 = — H = —gt 2 / 2, i.e. le temps de chute t = yj2H/g. Mais, d’autre part, 
Y — 0 = guj cos Xt 3 / 3 = (geo/3) cos A [2 H/g} 3 ^ 2 . A.Af. On trouve Y = +6, 5 cm. □ 
9. Nous avons vu que Y 1 s 2 = 5 10 -9 . 

10. La deuxième composante de R(t) — R 0 est positive pour tout t > 0 puisque 0 < A < n/2 
dans cet hémisphère. 
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3.2.7 Exemple : le pendule de Foucault (1819-1868) 


L’expérience du pendule de Foucault met en évidence de manière spectaculaire 
les effets de la force de Coriolis — due à la rotation diurne de la terre — sur le 
mouvement d’un pendule sphérique, par exemple le pendule de 67 m exposé au 
Panthéon (Paris). 

“Vous êtes invités à venir voir tourner la Terre. 1 ’ 


Léon Foucault (1851) 


Pour mettre en œuvre l’expression (3.2.26) de l’accélération relative du pendule 
de Foucault, nous négligerons comme précédemment les termes du second ordre en 
la vitesse angulaire u de la terre ; nous supposerons, de plus, que le pendule effectue 
des oscillations de faible amplitude (petits mouvement) dans un plan Z = const. 

Si nous ne tenions pas compte des forces d’inertie, le mouvement du pendule, de 
longueur I dans le champ de pesanteur terrestre g = const. > 0, serait régi par les 
équations du mouvement d’un pendule oscillateur harmonique bidimensionnel, i.e. 


X 

Ÿ 

0 



Ces équations du mouvement doivent maintenant être modifiées en utilisant (3.2.26) 

pour prendre en compte la vitesse instantanée de la terre (3.2.29) que l’on exprimera 

plutôt en terme de la latitude À du pendule ; on a alors 

/ X \ ( X \ ( D x = —a; cos À \ / X \ 

fU-f (y -2 U = o x (y . 

\ 0 / \ 0 J \ ^ sin A J \ 0 / 

Il vient alors, en posant Uq = yfgjï pour la pulsation propre du pendule, 


X-2ÇI z Ÿ + lü 2 0 X = 0, 
Ÿ + 2 DzX + üülY = 0 , 


ou encore, en posant W = X + ïY , 


W + 2 iü z W + ujIW 


0 . 


(3.2.35) 

(3.2.36) 


r\j 


(3.2.37) 
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L’intégration de (3.2.38) est aisée : on obtient 

W(t) = e~ lüzt (A cos u 0 t + B s'm u 0 t) (3.2.38) 


où A,BeC sont des constantes d’intégration. 11 

Le plan d’oscillation du pendule de Foucault possède donc une vitesse angulaire 
constante fl = — Ttz = —uj sin À par rapport à la terre. Ce plan effectue ainsi une 
rotation de 360° pendant le temps 


2vr _ Tq 
|fi| | sin À| 


(3.2.39) 


où T® désigne le jour sidéral (la période de rotation de la terre par rapport aux 
étoiles fixes). 

A.Af. Le plan d’oscillation du pendule de Foucault du Musée des Arts et Métiers 
(Paris) a donc pour période T = 23 h 56' 4"/ sin(48°50') = 31 h 48'. 


Remarque 3.2.14. Noter que f2<0sz0<À< 90° et donc que la rotation du plan 
d’oscillation du pendule s’effectue, dans l’hémisphère nord, dans le sens des aiguilles 
d’une montre. A l’équateur, O = 0 ? le pendule de Foucault oscille dans un plan fixe. 


11. Signalons que la solution générale de (3.2.38) est, en fait, W(t) = 
e -iQ. z t (^4 costal t + Bsinwif) avec A, B G C et u)\ = yTup + H 2 Z ; mais uj \ = wq(1 + 0(P|)) = ujq. 
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Chapitre 4 


Mécanique du solide 

4.1 Dynamique des systèmes 

Rappelons, afin de préparer à l’étude des équations de la mécanique de solide, les 
résultats fondamentaux concernant la dynamique des systèmes de points matériels en 
interaction mutuelle et soumis à l’action de forces extérieures. Nous nous proposons 
ici de rappeler les théorèmes généraux de la mécanique. 

Considérons un système de N points matériels Mi , . . . , Mm de niasses 
évoluant dans l’espace euclidien E 3 sous l’effet de forces d’interaction mutuelles 
fi2, • • • , fjv— i,iv et de forces extérieures f | Jxt , . . . , f ^ . Nous supposerons que le référentiel 
dans lequel les forces sont ainsi exprimées est un référentiel inertiel “fixe” . 

4.1.1 Théorème général I 

La force à laquelle est soumise le point M* est alors 

N 
3 = 1 

pour tout i — 1, . . . , N. La force totale agissant sur le système 

N N 

f=E £ ‘ = E f r‘ ( 4 -!-2) 

i — 1 i = 1 
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est, de fait, égale à la somme des forces extérieures en vertu de la troisième loi de 
Newton ou “Loi de l’action et de la réaction” 

f,/ - 0, 0 = (4.1.3) 

qui implique immédiatement % = 0- 

On rappelle que si p* = mv, désigne l’impulsion du point Mi à un instant donné, 
alors l’impulsion totale du système p = YliLi p, vérifie l’équation différentielle 

(4.1.4) 

où f est 1a. force totale extérieure (4.1.2). L’équation (4.1.4) résulte directement de 
la seconde loi de Newton : rrqr = f;, pour tout % — 1, . . . , N et constitue le Théorème 
général I. 

4.1.2 Théorème général II 

Ayant choisi un référentiel inertiel, donc une origine O, nous pouvons définir le 
moment de 1a, force au point M. t par k, = r j x f \ on r* = OM, pour tout i — 1, . . . , N. 
De même, le moment angulaire du point matériel M* , d’impulsion p, sera défini par 
ti = r* x pi pour font i = 1 , ,N. Ces quantités physiques dépendent 

explicitement de l’origine, O (on dit aussi point de base ) choisie. 

Le moment total des forces appliquées au système 

N N 

k = E k ' = Z r * xf r t (4-1.5) 

i = 1 i = 1 

est alors le moment total des seules forces extérieures. Ceci résulte encore de la 
troisième loi de Newton évoquée plus haut. 

Si i = ti désigne maintenant le moment angulaire total du système par 
rapport au point O, à un instant donné, l’évolution temporelle de cette quantité est 
gouvernée par l’équation différentielle 

(4.1.6) 
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où k est le moment total (4.1.5) des forces extérieures. L’équation (4.1.6) résulte 
aussi de la seconde loi de Newton et constitue le Théorème général II. 


Proposition 4.1.1. Lors d’un changement d’origine b = 00' 7 le moment angulaire 
se transforme comme suit 

£ = £' + b x p (4-1.7) 

où £ (resp. £' ) désigne le moment angulaire par rapport à O (resp. O' ) et p l’impul- 
sion totale du système. 

Démonstration. On trouve aisément 

N N 

£ = Y^ om < x p i = ^2 (°° + ° Mi ) x Pi = 00/ x P + £ ' 

1=1 1=1 

et donc le résultat attendu. □ 

Exercice 4.1.2. Prouver l’on a, de même, la loi de transformation suivante pour 
le moment total des forces 

k = k' + bxf. (4.1.8) 

Définition 4.1.3. On appelle torseur toute fonction M i— >• (L.P) de l’espace eu- 
clidien E 3 à valeurs dans M 3 x M 3 se transformant comme suit 

(L'.P') = (L + P x MM', P) (4.1.9) 


sous une translation MM'. 

Le couple (£, p) est appelé torseur cinétique et (k, f) représente quant à lui le 
torseur dynamique du système de points matériels considéré. 

Exercice 4.1.4. On appelle couple tout torseur dynamique (k, f) tel que f = 0; 
montrer que le moment des forces, k, est indépendant du point de base. 
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4.2 Configurations solides 

Un solide est, par définition, un ensemble de points matériels “indéformable”, 
c’est-à-dire dont tous les points restent à distance fixe les uns des autres au cours 
du temps. 

Définition 4.2.1. Un système de points matériels M 1; . . . , M N est appelé solide si 

||M,Mj|| = const. 

à chaque instant et pour tous i,j — 1, . . . , N. 

Nous étendons immédiatement cette définition au cas d’une distribution continue 
de masse (par exemple une toupie). 

Définition 4.2.2. Un système de points matériels S C E 3 constitue un solide si 

||MN|| = const. 

à chaque instant et pour tous M, N G S. 

4.2.1 Espace de configuration 

Comment maintenant fixer la configuration d’un solide ? 

Supposons donné un référentiel “fixe” euclidien 7 Z = (O, (e 1 e 2 e 3 )). Pour déter- 
miner la configuration d’un solide S par rapport à ce référentiel, il nous faut fixer 
trois points différents du solide, disons une origine O' G S et deux autres points 
M, N G S non tous trois alignés. Alors O'M et O'N sont indépendants. Soit e) 
la direction de O'M et e 2 celle de O'M x O'N ; la matrice A = (e) e 2 e 3 ) avec 
e 3 = e( x e 2 est, grâce à (3.1.10) une matrice de rotation, A G SO(3), complètement 
déterminée par 0',M,N qui fixent la configuration du solide. 

Nous venons de prouver qu’une configuration d’un solide S est déterminée par 
une matrice de rotation A = (e) e 2 e 3 ) définissant un repère orthonormé direct en 
un point O' lié au solide et par un vecteur b = OO' donnant la position du point O' 
par rapport à l’origine, O, du référentiel fixe. 1 


1. Une configuration d’un solide est donc, en définitive, repère euclidien TU = (O', (e( e 2 eg)). 
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Proposition 4.2.3. L’espace de configuration d’un solide est SO(3) x M 3 . 

Nous avons vu au chapitre précédent que cet espace n’est autre que le groupe 
(Spécial) euclidien, SE(3). 

4.2.2 Champ de vitesse dans les solides 

Nous pouvons tirer avantage des résultats obtenus dans le chapitre concernant 
les changements de référentiels non inertiels. 

Le solide S est maintenant animé d’un certain mouvement et sa configuration à 
l’instant t est ainsi déterminée par un couple (4(f),b(f)) G SE(3). 

Nous avons vu (cf. (3.2.16)) que la vitesse absolue, v, d’un point M G S du 
solide est donnée par v = v re i + b + u> x r rel où v re i = 0 est sa vitesse relative (nulle 
puisque le point M est fixe par rapport au solide), b = OO' et r re i = (PM ; ici uj 
représente le vecteur instantané de rotation du (référentiel lié au) solide par rapport 
au référentiel fixe. On résume par la formule donnant la vitesse d’un point M du 
solide, à savoir vm = vo> + uj x O'M. Les point M et O' étant arbitraires, nous 
avons prouvé la 

Proposition 4.2.4. La vitesse Vm d’un point M d’un solide et reliée à la vitesse vn 
d ’un autre point N par la relation suivante 

(4.2.1) 

où uj est le vecteur instantané de rotation du solide par rapport au référentiel fixe. 

Remarque 4.2.5. Le vecteur instantané de rotation uj ne dépend pas du choix 
d’une origine O' , il n’est défini que par la matrice de rotation A, cf. (3.2.4) ; l a ùù 
de transformation (4-2.1) montre alors que le couple (v, a ;) est un torseur au sens 
de la Définition 4-1-3 : nous l’appellerons torseur cinématique. 

Exercice 4.2.6. Prouver la propriété d’équiprojectivité de la vitesse 

(v M , MN) = (vzv, MN) 


vm = vm + MN x uj 


valable pour tous M, N G S . 
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4.3 Cinétique des solides 

La notion qui remplace, dans le cas continu, la distribution de masse (777^=1,...^ 
est celle de densité de masse, fonction continue g : S — > R + caractéristique de la 
composition physique du solide. 

Définition 4.3.1. On appelle masse du solide (S, g) la quantité 

M — J g(r) dV (r) (4.3.1) 

où dV( r) = dxdydz désigne l’élément de volume canonique de M 3 . 

Remarque 4.3.2. Nous ferons l’hypothèse qu’un solide S est une partie compacte 
et orientée de l’espace E 3 , donc de volume fini 


V = / dV( r) < +oo. 

J s 

4.3.1 Centre d’inertie 

Introduisons la notion importante de centre de masse d’un système (discret ou 
continu) de points matériels. 


Définition 4.3.3. On appelle barycentre R d’un ensemble (M.,, de 
points matériels, associé à une origine O, le vecteur 


R 


1 

M 


N 

Y m d 

i = 1 


(4.3.2) 


où Yi = OM; désigne le rayon vecteur du point M. , relativement àO et M = J^iLi m i 
la masse totale du système. 

Le barycentre d’un solide ( S , g) est, de même, défini par 



avec la définition (4-3-1) de la masse M du solide. 
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Définition 4.3.4. On appelle centre de niasse (ou encore centre d’inertie) 

d ’ un ensemble de points matériels le point G défini par 


OG = R (4.3.4) 

où R désigne le bary centre du système relativement au point O. 

Remarquons que le centre de masse G est, contrairement au barycentre, indépendant 
du choix d’une origine O. 

Exercice 4.3.5. Prouver, par exemple dans le cas discret, que le centre de masse G 
est uniquement défini par la relation 

N 

m, GMj = 0 (4.3.5) 

i = 1 

qui en constitue une définition alternative. 

Un corollaire important du Théorème général I (4.1.4) et de la définition (4.3.2), 
(4.3.3) du barycentre est donné par la 

Proposition 4.3.6. Le mouvement R(f) du barycentre d’un système de masse to- 
tale M est régi par l’équation différentielle suivante 

MR = f (4.3.6) 

où f désigne la force totale extérieure. 

Le mouvement du barycentre correspond au mouvement d’un point matériel de 
masse, la masse totale, soumis à une force représentée par la force totale extérieure. 

Exercice 4.3.7. (i) Déterminer la position du centre de masse G d’un demi-disque 
homogène de rayon R. (ii) Même question pour un hémisphère plein de rayon R. 

Démonstration, (i) Considérons le demi-disque S défini par x 2 + y 2 < R 2 et 0 < y ; 
si M désigne la masse de S et g = M/QttR 2 ) sa densité on a clairement R = (0, Y) 
où Y = {1/M) J s gydxdy = {2/i tR 2 ) /[o_r] x [ 07I -] rsin.0 rdrdO, i.e Y = 4R/(3n). 

(ii) Si l’hémisphère plein de rayon R est défini par x 2 + y 2 + z 2 < R 2 et 0 < z, 
on trouve R = (0, 0, Z) où Z = | R. □ 
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Exercice 4.3.8. Déterminer la position du centre de masse G d’un cône homogène, 
plein, de hauteur h et de rayon R. 

Démonstration. Considérons le cône S, dont la pointe est l’origine O, défini en 
coordonnées cylindriques ( r,9,z ) par 0 < 6 < 2n, 0 < r < Rz/h et 0 < z < h. 
L’élément de volume en coordonnées cylindriques étant dV = rdrdOdz , le volume 
de S est V = J s dV = R 2 h. Le barycentre du cône est alors R = (0,0 , Z) où 
Z = (1/E) f s z dV = (1/E) Jq zdz f^ h rdr J^dd et finalement Z = | h. □ 

4.3.2 Opérateur d’inertie 

Considérons maintenant un solide mobile autour d’un point fixe O. En ayant 
recours à (3.2.16), et au fait que b = 0 et v re i = 0, on voit que la vitesse absolue 
d’un point r de ce solide est 

v = u> x r. (4.3.7) 

Supposons, pour l’instant, le solide formé de N points distincts. Le moment 
angulaire d’un point matériel r de masse m par rapport à O se calcule aisément ; on 
trouve £ = mr x v = mr x (w x r) = —mr x (r x a;), c’est-à-dire 

£ = —mj(r) 2 u>. (4.3.8) 

D’où le lemme suivant donnant le moment angulaire total du solide. 

Lemme 4.3.9. (i) Le moment angulaire total par rapport à un point fixe O d'un 
solide formé de N points matériels (r est donné par 

(4.3.9) 

où I est un opérateur linéaire appelé opérateur d’inertie du système; c’est un 
opérateur symétrique qui prend la forme 

N 

I = (4.3.10) 

i=l 

N 

= ]N| 2 1— ùfi (4.3.11) 

i=l 



dans le cas discret. 
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(il) Dans le cas d’un solide (S, g), le moment angulaire £ du solide par rapport 
au point fixe O est toujours donné par (4-3.9), et l’opérateur d’inertie par 

ï = - J e (r) jirfdVir) (4.3.12) 

dV{ r) (4.3.13) 

dans le cas continu. 


= / 0(r) ||r|| 2 1 — 


r r 


Démonstration, (ii) Il suffit d’appliquer (4.3.8) pour calculer le moment angulaire £, 
de chaque point r, : pour i — 1, . . . , N. Le moment angulaire total du système étant 
£ = YliLi^h l a formule (4.3.10) suit. L’expression (4.3.11) est déduite de (3.2.8). 
L’opérateur est trivialement symétrique, 1 = 1, puisque la dyade r i rj est symétrique. 

(ii) Le moment angulaire total du solide est £ = f s (g( r) r x v) dV( r) et l’expres- 
sion (4.3.12) est alors déduite de (4.3.8). □ 


Exercice 4.3.10. Prouver que l’opérateur d’inertie (4-3.11) prend, en coordonnées 
cartésiennes, la forme suivante 

( + EÜi m ÉVi + z f) ~ EÜi m i x iVi ~ Eii m i x Ui ^ 

I= -EÏi mxiVi +T l ?=i m i( x i +zf) — EÜi myiZi (4.3.14) 

V - EU] m XiZi - Eir rrii + Y^=i ™É X Ï + vf) ) 

Exercice 4.3.11. Prouver que l’opérateur d’inertie (4-3.13) prend, en coordonnées 
cartésiennes, la forme suivante 

( + fs Q(y 2 + z 2 )dV —J s gxydV — f s gxzdV \ 


/ = 


— f s g xy dV 
\ — J s g xz dV 


+ J s g(x 2 + z 2 )dV — f s g yz dV 
— f s g y z dV + f s g(x 2 + y 2 )dV J 


(4.3.15) 


Corollaire 4.3.12. Le moment angulaire L d’un solide (S, g) par rapport à un point 
fixe O, exprimé dans le référentiel mobile A(t ) est donné par 


£ = AL 


(4.3.16) 
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c’est-à-dire 

( 4 . 3 . 17 ) 

où Cl est le vecteur instantané de rotation exprimé dans le repère mobile, u ; = ATI, 
et T = A L IA l’opérateur d’inertie dans le repère mobile. 2 

Démonstration. Il suffit d’utiliser les formules de passage du repère fixe au repère 
mobile, formules rappelées dans l’énoncé. □ 


L = ICI 


Remarque 4.3.13. L’opérateur d’inertie est un opérateur symétrique, il est donc 
diagonalisable. Ses directions propres sont appelées directions principales ou axes 
principaux d’inertie du solide. En pratique, on ramène le calcul de l’opérateur 
d’inertie I à celui de ses valeurs propres Ji,/ 2 ,/ 3 . 


Exercice 4.3.14. (%) Calculer l’opérateur d’inertie I par rapport au centre d’une 
boule homogène de masse M et de rayon R. (ii) Même question pour une sphère 
homogène de masse M et de rayon R. (iii) Même question pour un disque homogène 
de masse M et de rayon R. 


Démonstration. Réponse : (i) Pour une boule B\ homogène et de masse M la densité 
est g = M/^tiRQ ; utilisant le fait trivial 


’&k 


x dxdydz — I y dxdydz = z dxdydz = - 


14t t R 5 




! r3 


3 5 


( 4 . 3 . 18 ) 


on trouve aisément / = diag(/i, R, I\) avec I\ = | MR 2 . 

(ii) Pour une sphère S\ homogène de même masse, g = M/(4irR 2 ) et puisque 




( 4 . 3 . 19 ) 


il vient finalement / = diag(/i, I\, I\) avec I\ = | MR 2 . 


2. Nous ne considérerons désormais comme objet fondamental que l’opérateur d’inertie I propre, 
c’est-à-dire relatif au solide. 
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(iii) Pour un disque plein D 2 R homogène et de niasse M dans le plan z — 0, on a 
g = M/ (n R 2 ) et 

f x 2 dS = [ y 2 dS=-'ïïR i (4.3.20) 

M Jd I 4 

d’où / = diag(/i, /i, / 3 ) avec R = | MR 2 et J 3 = \MR 2 . Nous remarquons que, pour 
ce solide plan, on a 

7 3 = R + h 

dans les axes principaux d’inertie. Ce fait est général comme le montre clairement 
la formule (4.3.15). □ 

Exercice 4.3.15. Déterminer V opérateur d’inertie I d’un cylindre homogène de 

masse M , de hauteur h et de rayon R par rapport à son centre d’inertie. 

Démonstration. On trouve / = diag(/i, R, / 3 ) avec I\ = \AL{K 2 + |/i 2 ) et, d’autre 
part, J 3 = \M R 2 . □ 

Exercice 4.3.16. Calculer l’opérateur d’inertie I d’un ellipsoïde plein, homogène, 

de masse M et de demi-axes ( a,b,c ). 

Démonstration. Le bord de l’ellipsoïde, E, a pour équation x 2 / a 2 + y 2 /b 2 + z 2 /c 2 = 1. 
Posons x' = xja, y' = y/c et z' = z/c de sorte que dans ces nouvelles coordonnées 
on a x' 2 + y' 2 + z' 2 = 1, i.e. l’équation d’une sphère de rayon R = 1. On a donc 
h — fs Q(y 2 + É 2 ') dxdydz = g f B 3 abc(b 2 y' 2 + c 2 z' 2 ) dx'dy'dz' . En utilisant le résultat 
(4.3.18) avec R = 1 pour les variables x', y ' , z' on trouve R = g a6c(47r/15)(6 2 + c 2 ). 
On sait, d’autre part, que g = M/(^nabc) ; ceci entraîne R — | M(b 2 + c 2 ). Nous 
obtenons enfin / = diag (| M(b 2 + c 2 ), | M (c 2 + a 2 ), | M ( a 2 + b 2 )) . □ 

Proposition 4.3.17. Les valeurs propres (R, 1-2, R ) de l’opérateur d’inertie vérifient 
les inégalités triangulaires 

R < R + R, 7 2 < J 3 + / 1 , ^3 < R + 7-2, 

dont certaines deviennent des égalités dans le cas de solides plans. 
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Démonstration. On voit (4.3.15) que R + / 2 = f s g(y 2 + 2 z 2 + x 2 ) dE, donc /i + / 2 = 

Il + 2 f s g z 2 dV > I 3 . On a égalité si z = 0, i.e. dans le cas d’un solide plan. □ 

Définition-Théorème 4.3.1. On appelle rotateur tout solide dont les points sont 
distribués sur une droite, par exemple la droite x = y = 0. L'opérateur d’inertie 
d’un rotateur est 

I = diag(/i,/i,0). (4.3.21) 

Démonstration. Trivial. □ 


Exercice 4.3.18. Donner l’opérateur d’inertie Ig par rapport au centre d’inertie G 
en fonction de l’opérateur d’inertie Io par rapport au point fixe O en prouvant que 

I G = Io + Mj{ R) 2 (4.3.22) 

où R = OG désigne le barycentre et M la masse totale. 


Démonstration. Effectuons le calcul dans le cas discret. On a OM, = OG + GM, 

pour tout i = 1 , ,N de sorte que j(OM,) 2 = j(OG) 2 + j(GM,;)j(OG) + 

j(OG)j(GMi)+j{GMi) 2 . La définition (4.3.10) de Io donne alors 

N 

Io = -YmdjOMf ) 2 

i = 1 
N 

~Y m i [i(OG) 2 + j (GM 1 )j(OG)+i(OG)i(GM î )+i(GM ! ) 2 ] 

i = 1 

= —M j(OG) 2 + Ig 


grâce à la définition (4.3.5) du centre de masse, et en appelant M la masse totale. □ 
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On peut maintenant déterminer l’énergie cinétique du solide (par rapport au 
repère fixe) entrant dans la définition du lagrangien. Nous exprimerons cette quantité 
à la fois dans le repère fixe et dans le repère mobile lié au point fixe, O. 


Théorème 4.3.19. L’énergie cinétique du solide mobile autour d’un point fixe, O, 
est donnée par la forme quadratique suivante en la vitesse angulaire 


T = /«) = !(«,€), 

= t(n,in) = i(S2,L). 


(4.3.23) 

(4.3.24) 


Démonstration. On sait que T = è YliLi m i|| v i|| 2 dans le cas discret, par exemple. 
Grâce à (4.3.7), on trouve 

N 


T = X r< l 


i= 1 
N 


mj(u), Tj x (lo x r *)) 


i = 1 
N 




i = 1 


-(a;, 

2 V ’ 


' N 

S 

. i=l 




uj) 


-(u;, Iu>) 


avec la définition (4.3.10) de l’opérateur d’inertie. Le reste de la preuve découle de la 
relation (4.3.9) entre moment angulaire et vitesse angulaire. L’expression (4.3.24) de 
l’énergie cinétique en terme des quantités précédentes exprimées dans le référentiel 
mobile lié au solide suit du fait que c o = AD et I = AI A -1 , avec A G SO(3). □ 


L’énergie cinétique d’un solide évoluant dans l’espace sans point fixe, se calcule 
maintenant aisément ; elle est donnée par le 



62 


CHAPITRE 4. MÉCANIQUE DU SOLIDE 


Corollaire 4.3.20. L’énergie cinétique d’un solide est somme de l’énergie cinétique 
de translation du centre d’inertie. G, et de l’énergie cinétique de rotation autour du 
point G, i.e. 

(4.3.25) 

où M est la masse totale, v G la vitesse du centre d’inertie G du solide, / G est son 
opérateur d’inertie relativement à G et co sa vitesse angulaire. 

Démonstration. Travaillons, comme précédemment, dans le cas d’un solide formé 
de N points. L’énergie cinétique est T = | m *ll v *l| 2 °ù> (4-2.1), la vitesse 

du point Mi est q = + w x GM. ; pour tout i = 1, , N. On trouve alors 

aisément T — | Y2iLi m i (II v g || 2 + 2 (v G ,u; x GM,) + ||w x GM, : || 2 ) ou encore, en 
développant, T = \ M ||v G || 2 + (v G x u>, ^ =1 m GM, : ) + \ m i 11^ x GM ill 2 > 
en désignant par AI la masse totale. Le deuxième terme s’annule en vertu de (4.3.5) 
et le troisième donné par (4.3.23) avec O = G. □ 

4.3.4 Dynamique du solide 

Donnons, dans ce chapitre, les équations du mouvement d’un solide en présence 
de forces extérieures. Nous supposerons, pour simplifier, ce solide mobile autour d’un 
point fixe, O. Cette restriction, peu fondamentale, permet la description complète 
d’un certain nombre de systèmes mécaniques de première importance comme les 
toupies, les gyrocompas, les gyroscopes dédiés à la stabilisation des véhicules, etc. 
On déduit du Théorème général II, à savoir de l’équation (4.1.6), le 

Théorème 4.3.21. Soit k le moment des forces extérieures appliquées à un solide 
relativement à un point O ; l'évolution temporelle du moment angulaire £ de ce solide 
est gouvernée par l’équation différentielle 


T = ^ M || v g || 2 + ^(uj,I g uj) 


= k. 


d£ 

dt 


(4.3.26) 
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Compte tenu de la relation (4.3.9) entre moment angulaire £ et vitesse an- 
gulaire u> du solide on obtient dans le cas très particulier d’un solide à symétrie 
sphérique où / = diag(/i, R, IQ le 

Corollaire 4.3.22. En l’absence de forces extérieures, le moment et la vitesse an- 
gulaires d ’un solide à symétrie sphérique sont des constantes du mouvement , 3 

£ = const. & U) = const. (4.3.27) 

Démonstration. On déduit de (4.3.26) et de k = 0 que £ = const. Mais £ = I iw, 
donc ta = £/ 1\ = const. □ 

Remarque 4.3.23. Attention £ et co ne sont des vecteurs parallèles que dans le cas 
des boules et des sphères (solides SO(3 ) -invariants). 

Illustrons le Théorème 4.3.21 par un exemple mettant en évidence le mouvement 
de précession d’une toupie avec point fixe dans le champ de pesanteur terrestre. (La 
toupie “tombe” — comme tout corps dans le champ de gravitation — mais sa chute 
“libre” est de nature plus subtile et complexe que celle d’un simple point matériel.) 

Exemple 4.3.24. Considérons une toupie symétrique, R = R, de masse M , avec 
point fixe O. Cette toupie est plongée dans le champ de pesanteur terrestre g = const. 
On désigne par R = OG son barycentre, par R = ||R|| la distance du point O au 
centre de masse G et par u = R/i? la direction du barycentre. Dans l’approximation 
gyroscopique où le spin est dominant — on appelle “spin” ou vitesse angulaire propre 
le vecteur u sp i n = cu sp i n u de rotation instantanée de la toupie autour de son axe de 
symétrie. Dans cette approximation on a u; sp i n = const. et 

£ = lu spin (4.3.28) 

où I — /.s est le moment d’inertie par rapport à l’axe de symétrie de la toupie. 


3. Une toupie libre à symétrie sphérique tourne uniformément autour d’une direction fixe ! 
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On obtient, grâce à (4-3.26) et au fait que k = R x Mg = MRu x g, l’équation 
différentielle £ = Iû > s pin = Iu sp - m \i = —MR gxu, conduisant à l’équation différentielle 
suivante gouvernant le mouvement de l’axe de la toupie 


du 

dt 


— CiJprec X U 


MR 


ou 


OJ prec 


lu. 


(4.3.29) 


spin 


La direction de la toupie évolue donc au cours du temps selon la loi u (t) = 4(t)u 0 
où A(t) est une matrice de rotation autour de la verticale (la direction de g) et de 
vitesse angulaire constante de précession 

MRg 


^prec 


Ilü, 


spin 


où g = ||g|| est l’intensité de l’accélération de la pesanteur.' 


4.3.5 Lois de la statique 

Un solide est en équilibre statique si chacun de ses points a une position fixe, 
donc une vitesse nulle ; l’impulsion totale et le moment angulaire total du solide 
sont donc nuis, i.e. p = 0 et £ = 0. Ceci signifie, cf. les Théorèmes généraux (4.1.4) 
et (4.1.6), que la somme des forces extérieures agissant sur le solide et le moment, 
par rapport à un point O, de ces forces sont nécessairement tous deux nuis. 

Définition 4.3.25. Soient f la somme des forces extérieures et k le moment des 
forces extérieures relativement à un point O agissant sur un solide. Les conditions 
d’équilibre de ce solide sont 

(4.3.30) 

Remarque 4.3.26. Les conditions d’équilibre (4-3.30) sont, en fait, indépendantes 
du choix du point de base O. 

Quand il y a contact entre solides, les conditions d’équilibre doivent s’appliquer 
à chaque solide. Le contact entre deux solides en équilibre est garanti par la présence 


f = 0 & k = 0 


4. On a bien À(t)A{t) 1 = j(u prec ). 
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de forces de contact appelées encore forces de réaction entre solides. Les solides 
pouvant éventuellement être reliés entre eux par des fils rigides, il conviendra aussi 
de tenir compte des tensions des fils dans le bilan des forces de liaison entre solides. 
Les forces extérieures auxquelles est soumis chaque solide individuellement incluent 
donc nécessairement les forces de contact, de liaison, de frottement. 

Exercice 4.3.27. On considère une échelle double BiAB 2 dont les montants ABi et 
AB 2 , de masse négligeable, sont articulés autour d'une charnière A. Les montants 
sont reliés par un fil B\B 2 inextensible. Une personne de poids P se trouve au 
milieu, G, du montant ABi de l’échelle dont les pieds B 1 et B 2 sont en contact sans 
frottement avec le sol horizontal. Le système étant en équilibre, déterminer (i) la 
réaction Raj exercée en A par la tige AB 2 sur la tige AB\, (ii) les réactions du 
sol R B i, Rb 2 et (iii) la tension T b x du fil en fonction de P et de l'angle a entre les 
montants de l’échelle et l’horizontale. 

Démonstration. Traitons les deux solides AB\ et AB 2 séparément. Rappelons que 
la Loi de 1 action et de la reaction (4.1.3) implique R<Ai TR.a 2 = 0 cl T Sl +Tb 2 — 0. 
Les conditions d’équilibre (4.3.30) s’écrivent maintenant 


Rni + Ri?i + + P — 0 

(4.3.31) 

BiA x Ra, + BiG x P = 0 

(4.3.32) 

Ra 2 + Rb 2 + T b 2 = 0 

(4.3.33) 

B 2 A x Ra 2 = 0 

(4.3.34) 


On obtient aisément, à partir de (4.3.34), AB 2 x Ra 2 = 0, i.e. Ra 2 //AB 2 . 

Prenons O = Bi comme origine, e x comme la direction de B ] B 2 et e y comme 
celle de —P. On a 
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où Ra, Rb 1 , Rb 2 , T sont des constantes et P = ||P|| le poids de la personne. Il vient, 
grâce à (4.3.31), 

f — .R 4 cos a: + T = 0 

(4.3.35) 


— R a cos a + T =0 

+Ra sin a + Rb l — P = 0 


et, grâce à (4.3.32), 

2LPa sin ck cos a — -LP cos a: = 0 (4.3.36) 

en désignant provisoirement par L la longueur des montants de l’échelle. Enfin (4.3.33) 


(4.3.37) 


entraîne 

+Ra cos a — T =0 
— Ra sin a + Rb 2 = 0 
On déduit de (4.3.35) que T = Ra cos a & Rb 1 = P — PAsina et de (4.3.37) que 
Rb 2 = Pa sin a. Enfin, (4.3.36) entraîne 

P „ 3 P „ P _ P 


Ra ~ 


4 sin a ’ 

ce qui constitue le résultat attendu 


/, = 3P = P 

PBi ^ i aX b 2 ^ , 


T = 


4 tan a ’ 


□ 


4.4 Equations d’Euler & mouvements de Poinsot 

Nous étudions dans ce chapitre le mouvement d’un solide libre autour d’un 
point O, c’est-à-dire mobile autour d’un point fixe en l’absence de forces extérieures. 

4.4.1 Equations d’Euler 

Le théorème (4.3.26) nous apprend alors que le moment angulaire du solide par 
rapport à O et dans le repère “fixe” est l = const. 

Théorème 4.4.1. L ’ évolution temporelle du moment angulaire L d'un solide libre, 
par rapport à un point fixe O, est gouvernée, dans un repère mobile lié au solide, 
par 5 

(4.4.1) 

où le vecteur Cl représente la vitesse angulaire du solide. 


ClL T O 

— = L x 12 

dt 


5. Ce théorème est dû à L. Euler, “Theoria motus corporum solidorum” (1765). 
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Démonstration. Puisque t = AL = const. où, cf. (4.3.16), A{t) G SO(3) représente 
la configuration du solide à l’instant t, on a I = AL + AL = 0. Il s’ensuit que 
ÀA _1 AL + AL = j(oo)AL + AL = 0, grâce à (3.2.4). Mais (3.2.10) conduit, si 
oj = Ail, à 7lj(f2)7l _1 ÆL+7lL = 0, i.e. à 7l(j(f2)L+L) = 0 qui achève la preuve. □ 


Exprimons maintenant les équations d’Euler (4.4.1) dans un référentiel propre du 
solide, c’est-à-dire dans un référentiel lié au solide dans lequel l’opérateur d’inertie 
soit diagonal, 


h 


1 = 


(4.4.2) 


où Ii, 1-2 et J 3 représentent les valeurs propres (nécessairement positives) de J. Mais 
le moment angulaire L est, rappclons-le, relié à la vitesse angulaire il par L = I il. 
cf. (4.3.17), i.e. 

/ Li \ ( hih \ 

(4.4.3) 



Les équations (4.4.1) s’écrivent alors dans ce repère 


/ i iii i \ / i\ii\ \ / iii \ 

I I2ÙI2 J = j ^2^2 1 X I H2 J 

V IAh J \ J \il 3 J 


ou encore 

( /ifli = (I 2 — I 3 ) f^2^3 

< I 2 ÙI 2 = (h — h) O 3 O 1 

\ I3&3 — (h ~ I2) hli ^2 


(4.4.4) 


Le système d’équations non linéaire (4.4.4) constitue les équations d’Euler régissant 
l’évolution temporelle de la vitesse angulaire d’un solide libre autour d’un point fixe. 


Remarque 4.4.2. (%) L’intégration de ce système d’équations différentielles est 
délicate, dans le cas général d’une toupie asymétrique où ii,/ 2 et I 3 sont tous 
différents ( elle met en jeu les fonctions elliptiques), (ii) Dans le cas particulier 
d’une toupie à symétrie sphérique, I\ — I 2 — I 3 , on retrouve, par contre, le résultat 
élémentaire du Corollaire 4-3.22. 
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4.4.2 Exemple : la toupie symétrique 

Si une toupie possède un axe de symétrie privilégié, disons de direction e ' 3 , 6 on 
est en présence d’un système à symétrie cylindrique. L’opérateur d’inertie de 
cette toupie prend la forme (4.4.2) avec 

h = hÉ h • (4.4.5) 

Dans ces conditions, il vient trivialement grâce à (4.4.4) 


D 3 = const. (4.4.6) 

De même, on trouve D 3 = ((/i — / 3 )/ii)D 3 D 2 et D 2 = ((h — /i)//i)D 3 Di, i.e. 
Di = —u)fl -2 et D 2 = ou encore, en posant D = D 3 + iD 2 , 


Ù = iuD avec 

La solution générale de (4.4.7), à savoir 



(4.4.7) 


ü(t) = D 0 e iuj[t ~ to] 


(4.4.8) 


avec D 0 = const. > 0 et f 0 G M met en évidence le fait que la projection D de la 
vitesse angulaire sur la plan orthogonal à l’axe de la toupie tourne avec une vitesse 
angulaire constante u donnée en (4.4.7). Grâce à (4.4.5) on conclut que la vitesse 
angulaire il et le moment angulaire L ont un mouvement de précession autour 
de l’axe de la toupie ; ils tournent autour de cet axe avec la même vitesse angulaire u. 


4.4.3 Mouvements de Poinsot 

Etudions maintenant le mouvement du moment angulaire d’un solide générique, 
par exemple d’une toupie asymétrique dont l’opérateur d’inertie I a des valeurs 
propres ii,/ 2 ,/ 3 différentes, par exemple telles que 

I\ > I 2 > / 3 . (4.4.9) 

Nous avons immédiatement un résultat général donné par le 


6. Rappelons que la configuration du solide est donnée par la matrice de rotation A = (ej e' 2 e' 3 )- 
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Théorème 4.4.3. Les équations d’Euler (4-4- V possèdent deux intégrales premières 

T 2 r 2 T 2 

2 E — (L, Q) = VI + LA + TA = const . (4.4.10) 

h h h 

L 2 — (L, L) = L\ + L\ + L\ = const. (4.4.11) 

Démonstration. L’énergie cinétique du solide (4.3.24) donnant ici l’énergie totale 
E = |(0, TEL) = |(0, L), on obtient d.E/dt = (ft.l dfl/dt) = (Q,dL/dt) et 

dE/dt = (O.L x O) = 0 grâce aux équations d’Euler (4.4.1). On a, de même, 

d{L 2 )/dt = 2(L. dL/dt) = 2(L. L x O) = 0 grâce à (4.4.1). □ 

Exercice 4.4.4. Prouver (4-4-10) et (4-4-H) en utilisant la forme (4-4-4) des 
équations d’Euler. 

Nous voyons donc que le mouvement du moment angulaire L s’effectue sur des 
courbes, intersection d’un ellipsoïde de demi-axes \f2EÎ[ ï \J2 EI 2i \j2EI :i et d’une 
sphère de rayon L 2 . Les caractéristiques de ses surfaces sont, bien sûr, définies par 
les conditions initiales que l’on impose au système. Ces courbes sont clairement des 
courbes fermées, impliquant que le mouvement du solide est périodique. 

Six points jouent un rôle particulier; ce sont les sommets de l’ellipsoïde corres- 
pondant à des valeurs de l’énergie, E, et du moment angulaire, L , telles que les 
courbes se réduisent à des points (par exemple L = \J2 EIQ. Ces sommets corres- 
pondent donc à des valeurs constantes du moment angulaire L qui est donc, comme 
la vitesse angulaire fi, une constante du mouvement : il existe, pour un solide libre 
autour d’un point fixe O six mouvement de rotations stationnaires autour des 
trois axes d’inertie du système. 

Pour d’autres relations entre E et L il existe des courbes fermées correspondant 
à des trajectoires (périodiques) de L. Dans le cas (4.4.9), on montre qualitativement 
en étudiant ces courbes aux voisinages des sommets de l’ellipsoïde que les rotations 
stationnaires autour des directions extrêmes ej et e' 3 sont stables 1 alors qu’elles sont 

7. Les courbes fermées restent dans un voisinage des sommets qu’elles entourent quand on 
diminue par exemple la valeur du moment angulaire L. 



70 


CHAPITRE 4. MÉCANIQUE DU SOLIDE 


instables 8 autour de la direction médiane e' 2 . Voir la Figure suivante. 9 



Figure 4.1 - Toupie de Poinsot 


Afin d’étudier maintenant le mouvement du solide relativement au repère fixe, 
introduisons un objet géométrique utile, l’ellipsoïde d’inertie 


£ = {X G R 3 |(X,XX) = 1} 


(4.4.12) 


relativement au repère mobile. Rappelons, cf. (4.4.10), que 2 E = (f2,Xfî), donc 


X 


n 

V2Ë' 


(4.4.13) 


Si A(t) = (e( e' 2 e' 3 ) G SO(3) repère la position du solide autour du point O dans 
le référentiel fixe au temps t, on désignera par £ t = A{t)£ la position de l’ellipsoïde 
d’inertie au temps t par rapport au repère fixe. 


Théorème 4.4.5 (Poinsot). L’ellipsoïde d’inertie Et roule sans glisser sur un plan 
fixe perpendiculaire au moment angulaire constant L. 

Démonstration. Le vecteur N = grad(X,XX) = 2 XX est automatiquement ortho- 
gonal à la surface (4.4.12) au point A" G £. Nous allons montrer que la norma- 
le n = A(f)N à £ t au point x = A(f)X est colinéaire au vecteur (constant) £. 
En effet, n = 2 A(t)XX = 2A(t)T Ut/ y/2E grâce à (4.4.13). Mais, cf. (4.3.17), le 

8. Les courbes s’éloignent immédiatement des sommets médians après une modification infini- 
tésimale du moment angulaire L. 

9. http : / / www . cds . caltech . edu/“marsden/books/Mechanics_and_Symmetry . html 
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moment angulaire est donné par L = Z Tl relativement au solide et par £ = A{t)L 
dans le repère fixe, cf. (4.3.16) ; il s’ensuit que n = 2£j\/TË est bien parallèle au 
moment angulaire £. Le plan tangent T^£ t en x G £ t est donc perpendiculaire au 
vecteur constant £. Montrons que ce plan reste fixe au cours du temps. On a (£, x) = 
(L, X) = (L. fl)/\/2E, grâce à (4.4.13) ; mais (4.4.10) entraîne alors (£,x) = \J2E = 
const. La projection du vecteur x(t) G £ t sur la direction du vecteur fixe £ reste 
constante : le plan T x £ t est indépendant du temps t. Comme x = uxx = 0 puisque 
x est parallèle au vecteur instantané de rotation, u>, la vitesse du point de contact 
x G £ t avec le plan fixe T x £ t est nulle et l’ellipsoïde d’inertie roule sans glisser. □ 

Exercice 4.4.6. (i) A quelle condition f2 0 (t) = (0,0,f2o(£)) est-elle une solution des 
équations d’Euler (4-4-4) pour une toupie asymétrique ? (ii) Linéariser les équations 
d’Euler au voisinage de cette solution, (iii) En déduire l’expression a-p-prochée de 
l’évolution temporelle du moment angulaire L(t) = L 2 (t), L 3 (t)). 

Démonstration, (i) Supposons que l’on ait (4.4.9) avec I 3 > 0. Alors, nécessairement 
Oo = 0 et donc O 0 = const. (rotation stationnaire autour du troisième axe). 

(ii) Posons f î(t) = f2 0 (t) + eZ (t) avec e < 1 ; on obtient Z\ = Z- 2 £l 0 (l 2 — h) /h 
et Z -2 = ZiQ q (I 3 — /i)// 2 avec Z 3 = 0 modulo des termes en 0(e 2 ). 

(iii) Revenant au moment angulaire L = (RLli, I- 2 fl 2 , 1 3 0 3 ), nous avons donc 



Cherchant la solution générale de ce système avec l’Ansatz général : L\ = Aie lult et 
L 2 = A 2 e lut avec wGlet Ai, A 2 G C, on trouve 



avec A 2 = iA\\J (1 — I 3 /Ii)/(1 — / 3 // 2 ). Posons A\ = Lyjl — I 3 // 2 où L est une 
constante d’intégration que l’on peut choisir réelle ; on trouve alors 

Li(t) = lJ 1 — -J- cosc ot, L 2 (t) 

Le moment angulaire L(t) décrit donc une petite ellipse autour du vecteur e' 3 . □ 


^ —L 


1 — y - sin eut, 
Ii 
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4.5 Toupie de Lagrange 

Nous abordons enfin, dans ce chapitre, l’étude qualitative des mouvements de la 
toupie de Lagrange, c’est-à dire d’une toupie, à symétrie axiale (ou cylindrique), 
mobile autour d’un point fixe dans le champ de pesanteur constant. 

La toupie (S, g) supposée de densité de masse constante, g = const., est donc 
mobile autour d’un point fixe, sa pointe O, dans le champ d’accélération de la 
pesanteur g = const., d’intensité g = ||g||. On désigne par G le centre d’inertie du 
solide, par R = OG son barycentre relativement à O et on pose R = ||OG||. 

4.5.1 Angles d’Euler 

Comment repérer une configuration de la toupie ? En d’autres termes, comment 
paramétrer la matrice A G SO(3) repérant la position de la toupie, à un instant 
donné, relativement à un repère euclidien fixe 1Z = (O, (ex e 2 e 3 )) ? 

La réponse à cette question est due à L. Euler qui a proposé de cartographier le 
groupe SO(3) des rotations euclidiennes à l’aide de trois angles appelés depuis lors 
angles d’Euler. 

Considérons les directions propres de la toupie confondues dans un premier temps 
avec la base orthonormée directe fixe (ex e 2 e 3 ). Pour atteindre une configuration 
arbitraire de la toupie procédons en trois étapes. 

1. Line rotation A 3 {<fi) d’angle (f> autour de e 3 transforme la base (e 3 e 2 e 3 ) en 
(e* e 2 e 3 ) = (ex e 2 e 3 )4 3 (0) en invariant e 3 = e 3 . 

2. Line rotation AfiO) d’angle 6 autour de e* transforme la base (e* e 2 e 3 ) en 
(eî* e” e 3 *) = (ej e* 2 e£)Ai(0) en invariant ej = ej*. 

3. Une rotation A 3 (x/j) d’angle t/j autour de e 3 * transforme la base (e 3 * e 2 * e 3 *) 
en (e'x e' 2 e 3 ) = (e” e 2 * e 3 *)4 3 ('0) en invariant e 3 * = e 3 . 

Nous avons enfin (e'x e 2 e 3 ) = (ex e 2 e 3 )A 3 ((p)Ai(9)A 3 (fi>) ; en identifiant comme 
convenu (ex e 2 e 3 ) à la base canonique de M 3 , l’expression de la base (e'x e 2 e 3 ) = A 
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attachée à la toupie prend la forme suivante 

A = A 3 (<f>)A 1 (e)A 3 ^) ( 4 . 5 . 1 ) 

Définition 4.5.1. L'angle 0 < (j) < 2n est appelé azimut et 0 < 9 < n angle de 
nutation, l’angle 0 < if < 2 tt étant quant à lui l’angle de rotation propre (ou 
spin) ; voir la Fig. f.2. 10 L’axe engendré par est appelé ligne des nœuds. 



Les trois angles d’Euler fournissent ainsi une “carte locale” 

]0, 27t[x]0, 7t[x]0, 27t[ — » SO(3) : (0, 9, ip) i— >• A(4>, 9, iff) 

du groupe des rotations. (Les rotations infinitésimales forment, cf. (3.2.4), un espace 
vectoriel réel de dimension 3 constitué des vecteurs instantanés de rotation u) G M 3 .) 
Nous avons donc la 

Proposition 4.5.2. Les matrices A G SO(3) sont de la forme générique 

A{4>, 9, if) = exp (0j(e 3 )) exp(0j( e i)) exp(^j(e 3 )). (4.5.2) 

Démonstration. La preuve découle du résultat (4.5.1) et du fait qu’une rotation 
d’angle a autour d’une direction fixe u est de la forme A = exp(o j(u)) où l’exponen- 
tielle d’une matrice Z est définie par la série 

y2 yn +°° yn 

exp(Z) = l + Z+ — + — + ...=^ — 

n=0 

dont on démontre qu’elle est convergente. □ 

10. http : //mathworld. wolfram. com/EulerAngles .html 
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Exercice 4.5.3. Vérifier que Von a bien 

/ 1 0 0 \ / cos — sin fi 0 \ 

Ai(6) — I 0 cos 9 — sin 9 J & A 3 (fi) = I sin fi cos fi 0 J 
y 0 sin 6 cos 9 J y 0 0 1 J 

4.5.2 Lagrangien de la toupie de Lagrange 


Déterminons le lagrangien d’un solide quelconque conduisant aux équations régis- 
sant son mouvement. Il nous faut donc déterminer l’énergie cinétique T et l’énergie 
potentielle V en ternie des angles d’Euler (fi, 9, fi) et des vitesses associées (fi, 9, fi). 
Nous spécialiserons ensuite ce lagrangien au cas de la toupie de Lagrange. 

L’expression (4.3.24) de l’énergie cinétique met en jeu la vitesse angulaire VI 
du solide exprimée dans le repère mobile 7 Z' = (O, (ej e' 2 e' 3 )) par rapport auquel 
l’opérateur d’inertie est diagonal, i.e. 

( h 

X= I 2 

\ h 

où les moments d’inertie /i,/ 2,/3 du solide sont a priori arbitraires. La toupie de 
Lagrange est, elle, caractérisée comme en (4.4.5) par. 

(4.5.3) 


h = h + h 


Lemme 4.5.4. La vitesse angulaire du solide exprimée dans IZ 1 est donnée par 

( fi sin 9 sin fi + 9 cos fi \ 


VI = 


fi sin 9 cos fi — 9 sin fi 
\ fi cos 9 + fi 


(4.5.4) 


/ 


Démonstration. On a A~ l À = A~ 1 (ÀA~ 1 )A = A~ 1 j(u>)A avec la définition (3.2.4) 
du vecteur instantané de rotation eu. De plus (3.2.10) entraîne A~ l À = j(A~ l uj), ou 
encore 


A~ l À = j(Vt) 


(4.5.5) 


grâce à (3.2.18). Cette formule donnera, pour une matrice A de la forme (4.5.1), 
l’expression (4.5.4) recherchée. 
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On trouve, en effet, À = À 3 ({f>)A 1 (9)A 3 (-ip)+A 3 ((p)À 1 (9)A 3 (ip)+A 3 ((p)A 1 (9)À 3 (ip) 
et donc, puisque A ~ 1 = A 3 {^>)~ 1 AiiO) -1 A 3 {(f))^ 1 , il vient 

A 1 À= + A 3 (ip) l Ai{6) 1 A 3 ((fi) 1 À 3 {(p) A 3 {9) A 3 {y>) 

+ A 3 (xJj) l Ai{6) 1 Ài{6) A 3 (ip) 

+ A 3 {^)~ l À 3 {^) 

= + A 3 {'ii)y l A l (ey 1 (t)j(e 3 )A l (e)A 3 {i )) 

+ A 3 (yy l è jy^Ayy 

+ Vote s) 

= + 4>j(A 3 (-'ijj)A 1 (-6)e 3 ) 

+ ^ j {A 3 {—^)ei) 

+ VO( e 3) 

= + (p j(sin6 siu'i/j e 3 + sin d cos ^62 + cos de 3 ) 

+ 6? j (g°s e i — sin. ^ e 2 ) 

+ Vo( e s), 

c’est-à-dire l’expression attendue. □ 

L’énergie cinétique du solide T = | (O, TEL) = \ (/iO^ + I2D-2 + J 3 f2§) est alors 

T— + — ((^sindsin^ + dcos'V ) 2 
2 

-^2 * * 9 

+ — (^sinÆcos^ — dsimjj) (4.5.6) 

+ ^-ÿœsd + îy 2 . 

L’énergie potentielle du solide (5, g) est V — f s ggzdxdydz = Mg(R., e 3 ), i.e. 

V = MgR cos 9. (4.5.7) 

Proposition 4.5.5. Le lagrangien de la toupie de Lagrange (4-5.3) est donné par 

(4.5.8) 
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Démonstration. Nous avons L = T — V où T est donnée par (4.5.6) et V par (4.5.7). 
Le cas particulier I\ = / 2 , cf. (4.5.3), conduit alors au lagrangien L recherché. □ 


4.5.3 Mouvements de la toupie de Lagrange 

Nous sommes maintenant en mesure de présenter une étude qualitative des 
mouvements de la toupie de Lagrange en tirant proût de l’existence de trois intégrales 
premières des équations du mouvement. 

Proposition 4.5.6. Les équations de Lagrange associées au lagrangien (4-5.8) pos- 
sèdent trois intégrales premières, à savoir 


Q = sin 2 6 + I3 cos 2 9) + I^ip cos 9 = const. (4.5.9) 

Z/ 3 = J 3 (0 cos 6 + Tp) — const. (4.5.10) 

H = ^-(0 2 sin 2 6 + 9 2 ) + y (<j>cos9 + ip) 2 + MgRcos9 — const. (4.5.11) 

Démonstration. Les variables 0 et ^ sont clairement des variables cycliques (1.3.17) 
car le lagrangien L n’en dépend pas. Les quantités £3 = dL/dcp et L 3 = dL/d'ip sont 
donc des constantes du mouvement. Remarquons que £3 = (£, e 3 ) est, en fait, la 
troisième composante du moment angulaire £ ; quant à L 3 = / 3 fl 3 , c’est la troisième 
composante du moment angulaire L exprimé dans la base mobile. Enfin H = T + V 
n’est autre que l’hamiltonien, constante du mouvement car, cf. (2.1.6), on a dH/dt = 
—dL/dt = 0. □ 


Nous obtenons immédiatement les équations différentielles suivantes pour ip{t) 
et <f>(t), résolubles par quadratures 11 une fois connue la fonction 9{t), à savoir 


ainsi que, grâce à (4.5.9), 


L 3 

ip = — (p cos 9 

h 


£3 — L 3 cos 6 


h sin 2 9 


(4.5.12) 


(4.5.13) 


11. On dit qu’un système d’équations différentielles est résoluble par quadratures si l’intégration 
de ce système se ramène à un simple calcul de primitives. 
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L’hamiltonien (4.5.11), compte tenu de (4.5.10) et (4.5.13), prend alors la forme 


h=!ï 

2 


I 3 — L 3 cos 9 
h sin 2 9 


sin 2 9 + 9'- 


+ — ( — J + AI g R co s 9 


qui conduit à la 


Proposition 4.5.7. L'angle de nutation 9(t ) de la toupie de Lagrange varie en 
fonction du temps t de telle sorte 

E= I ±9 2 + V eS {9) = const. (4.5.14) 

avec la définition suivante du potentiel effectif 

V eS (9) = — + MgRcos9. (4.5.15) 

21 1 sin 9 

Démonstration. Il suffit de développer l’expression précédente de H et de poser 
E = H — L§/(2/ 3 ) qui est, comme H et L 3 , une constante du mouvement. □ 


Nous allons utiliser l’intégrale première de l’énergie E = const. du système ainsi 
réduit à un degré de liberté pour déterminer qualitativement les différents mouve- 
ments possibles de l’axe de la toupie. 

Définissons les constantes suivantes 


_ I 3 r L 3 _ 2 E 

a — , b —— , ol — , 
D D 


fi 


2AIgR 

h 


> 0 


et posons 


u = cos 9. (4.5.16) 

Il vient alors ù = —sin 9 9 et donc 9 2 = ü 2 /(l — u 2 ) ; ainsi (4.5.14) s’écrit, grâce 
à (4.5.15), comme E = \l\Ct = \l\ [h 2 /(l — u 2 ) + (a — bu) 2 /{I — u 2 ) + fiu\. Nous 
obtenons finalement l’équation suivante 


u 2 = f{u) où f{u) = (a — fiu){l — u 2 ) — (a — bu) 2 (4.5.17) 


avec 


a — bu 
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ainsi que 


il, u(a — bu ) 

É = t b 


(4.5.19) 


7.3 1 - M 1 2 

Notons que le polynôme (4.5.17) du troisième degré f(u) = Pu 3 + . . . tend vers 
±cx) quand u — > ±oo. D’autre part, (4.5.16) imposant — 1 < u < 1, on a en fait 
/(±1) = —(a. =F b) 2 < 0, si a ^ ±b , aux bornes de l’intervalle permis pour u. Il 
est donc clair que le polynôme f{u) a, pour un mouvement réel de la toupie, trois 
racines réelles mi ,«2 et «3 telles que, génériquement, — 1 < u± < 112 < +1 < «3. Mais 
ù 2 = f{u) > 0 impose enfin u G [iq , m 2 ] • L’angle 0(f) que fait l’axe de la toupie avec 
la verticale varie donc périodiquement entre deux valeurs 9 1 et 62 (avec cos 9k = Uk 
pour k — 1,2) dépendant des conditions initiales, c’est le mouvement de nutation 
de la toupie 

62 < 6{t) < di. 


Enumérons maintenant, grâce à (4.5.18), les différentes lois horaires p(t) de l’azi- 
mut, c’est-à-dire les mouvements de précession. Posons 


, a 

U ~ b ~ L-i 


1. Si u' $:]«!, «2 [, alors p(t) 7^ 0; la fonction p(t ) est monotone et l’axe de la 
toupie effectue un mouvement de précession uniforme autour de la verticale. 


2. Si u' g]wi , U2 [, alors p(t) change de signe au cours du temps t; la direction de 
l’axe de la toupie effectue des aller-retours selon des courbes en feston sur la 
sphère unité. 


3. Si u' = U\ (resp. u' = U2), la vitesse 4>(t) s’annule périodiquement sur le 
parallèle 6 = 61 (resp. 6 — 62)', la direction de l’axe de la toupie décrit une 
courbe avec pointes (points de rebroussement) sur la sphère unité. 


Exercice 4.5.8. Il existe une solution particulière des équations de Lagrange avec 
axe de la toupie vertical, 6{t) = 0 ; cette solution, telle que !■$ = L 3 = /3D3 = const. 
avec D3 = 0 + ip est appelée toupie dormante. A quelle condition la position 
d'équilibre 9 = 0 est-elle une position d’équilibre stable ? 
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Démonstration. Développons le potentiel effectif (4.5.15) an voisinage de 6 = 0 en 
posant 9 = e9\ avec £ -C 1. Il vient aisément 



et 6 = 0 est une position d’équilibre stable si V^(0) > 0, c’est-à-dire si 


n2 ^ 4M g RI ! 

> 2 

1 3 


En présence de frottements, la vitesse angulaire fl 3 diminue et viole l’inégalité ci- 
dessus : la toupie se . . . réveille. □ 
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